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1. Introduction 



The present notes are an effort to explain in relatively non-technical terms re- 
cent results in many-body scattering and related topics. Thus, many results in the 
many-body setting should be understood as new results on the propagation of sin- 
gularities, here understood as lack of decay of wave functions at infinity, with much 
in common with real principal type propagation, i.e. wave phenomena. Motivated 
by this, I first briefly describe propagation of singularities for the wave equation. 
This is a remarkable relationship between geometric optics (the particle view of 
light) and the solutions of the wave equation (the wave view). 

Next, in Section I explain the geometry of many-body scattering, which in- 
cludes both that of the configuration space and phase space. This geometry is 
closely related to classical mechanics, playing the role of geometric optics, but even 
at this point quantum phenomena emerge. This leads to the analytic results, namely 
the propagation of singularities connecting classical and quantum mechanics. 

Much as for the wave equation, such a result has immediate applications, includ- 
ing the description of the scattering matrices and of the scattering phase. Slightly 
stronger versions can even lead to inverse results, a topic covered in the following 
section. 

After so explaining the results, in Sections [|^, I will try to at least give a flavor of 
how they are proved. This uses a many-body pseudo-differential algebra and posi- 
tive commutator estimates, so these are discussed. We remark that these techniques 
are closely related to the proofs of the propagation of singularities for the wave equa- 
tion, but there are significant differences as well, mostly arising from bound states 
of particles, which have no analogues for the wave equation. The pseudo-differential 
algebra itself is very interesting from the viewpoint of non-commutative geometry: 
there is a hierarchy of operator valued symbols at infinity. 

Asymptotic completeness was the main focus of work in many-body scattering 
for a long period. In Section [7], I briefly explain how it relates to the microlocal 
estimates. 

There is another area that is very closely related to many-body scattering, 
namely scattering on higher rank non-compact symmetric spaces. Here, in Sec- 
tion ||, we only discuss rank two, which corresponds to three-body scattering, since 
this is the only part that has been properly written up, but it is expected that very 
soon these results will extend to all higher rank spaces. 

I hope that the notes will make many of these results more accessible, the con- 
nections more transparent, and explain the motivation behind them. Many-body 



This work is partially supported by NSF grant #DMS-0201092, a Fellowship from the Alfred 
P. Sloan Foundation, and the Universite de Nantes, where these lectures were originally given. 

1 



2 



ANDRAS VASY 



scattering has a long history, and here I can only talk about the most recent devel- 
opments. An excellent overview of results known in the early 1990s can be found 
in Hiroshi Isozaki's lecture notes |3C|j . Indeed, in some sense, the current notes 
continue where J30| left off. I introduce a fully microlocal picture, motivated by 
the geometric approach of Richard Melrose p2|] , and emphasize the results these 
give, but the basic spectral and scattering results follow from a simpler 'partial' 
microlocalization, which is one of the subjects of |5Gj. 

The notes were originally prepared for a mini-course at the Universite de Nantes 
at the invitation of Professor Xue-Ping Wang, whose hospitality I gratefully ac- 
knowledge. The analytic continuation of the resolvent on symmetric spaces is a 
more recent development, but it was fueled by a discussion during the visit of Rafe 
Mazzeo, my collaborator, to Nantes. Over an espresso, Gilles Carron mentioned 
that the existence of the analytic continuation was not known, something that was 
hard to believe, but we immediately realized that our methods should yield such 
a continuation rather directly. I also thank Gunther Uhlmann for urging me to 
write up these notes: without him, they may never have been written up, and Rafe 
Mazzeo for a careful reading of the manuscript. 



2. Geometric optics and the wave equation 

According to the rules of geometric optics, light propagates in straight lines, and 
reflects/refracts from surfaces according to Snell's/Descartes' law. That is to say, 
considering light as a stream of billiard balls, the energy as well as the tangential 
component of the momentum (tangential to the surface hit) is conserved upon 
hitting the surface. 

But light satisfies the wave equation, i.e. if u — u{x, t) is the electromagnetic 
field on Vi x x M. t , Vl C R n , then Pu — where P is the wave operator c 2 A — Df, 
and a boundary condition also holds (say, Dirichlet), if 51 is not the whole space. 
(Here Dt = i and A = ^ ■ D\ is the positive Laplacian.) How are these two 
viewpoints related? 

One can phrase the connection in different ways. The most usual one in physics 
is that the billiard ball picture is accurate in the high frequency, i.e. low wave 
length, limit. That is to say, for high frequency light, geometric optics is accurate 
up to a 'small' error. A slightly different way of looking at this, which however does 
not involve approximations, is that the location of singularities of the solution of 
the wave equation is exactly predicted by geometric optics. Here singularities are 
understood as lack of smoothness, or possibly lack of analyticity. 

Indeed, it is convenient at this point to generalize the setting somewhat. So let 
(£l,g) be a Riemannian manifold with corners, P — c 2 A g — D 2 , c > 0. The speed 
of light, c, may be absorbed in the metric g, of course, we keep the notation in 
analogy with the usual wave equation. 

For simplicity of notation in this paragraph we assume that M z = fl x x R t is 
boundaryless; in general, the same definitions hold in the interior of M. Thus, we 
associate a homogeneous real function on T* M to P, namely its principal symbol: 
P = c2 l£lg — r2 > wnere we write C = (£, T ) as the dual variable of z = (x,t). Now 
T*M is a symplectic manifold with symplectic form u> = J2^Cj A dzj. Thus, p 
gives rise to a vector field H p , called the Hamilton vector field, by requiring that 
u>(V, Hp) — Vp for any vector field V on T*M. Hence H p is a smooth vector field 
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Figure 1. Projection of broken bicharacteristics to fi. When rays 
hit the boundary hypersurfaces, the tangential component of the 
momentum and the kinetic energy are conserved, but the normal 
component may change. At the corner, there is no tangential com- 
ponent (though there would be if the time variable were not pro- 
jected out), so the only constraint is the conservation of kinetic 
energy. 



on T*M explicitly given by 

dp d dp d 
p= dC ~d~z~ dz"d(' 

Note that p is constant along the integral curves of H p since taking V = H p , 
= u>(H p , Hp) — H p p. Null bicharacteristics are the integral curves of H p inside 
its characteristic set E = p~ 1 ({0}). Thus, if 7 : 1 — * T*M is a null bicharacteristic 
(here / is an interval), and z(s) = z(~f(s)), ((s) = (("f(s)), then these solve the 
ODE's % = §|, g = -ff . Hence, when M = O x R, Q c R", p = c 2 |£| 2 - t 2 
as above, we deduce that ^ and r are constant along the integral curves of H p , 
hence their projection to M consists of straight line segments. More generally, the 
projection of null-bicharacteristics to f2 are geodesies of g. 

There is an appropriate extension of this at boundary surfaces and even at cor- 
ners, called generalized broken bicharacteristics, see |36| , which I will not explain 
in full generality, though I remark that many-body scattering, discussed in the next 
section in detail, is rather similar. However, a somewhat typical example is that 
of broken bicharacteristics. These are piecewise bicharacteristics, i.e. there is a 
sequence Sj, j in a subset of integers, such that for each j, j\r s -,a j+ i) ^ s a nu ^ 
bicharacteristic in the sense described above, the projection z o 7 of 7 to M is con- 
tinuous, and j(sj+) — j(sj—) is conormal to the smallest dimensional boundary 
face containing z(j(sj)). Thus, the tangent vectors to zo<y\t a s \ and zo^\^ s ,_ ua .\ 
differ by a vector normal to the smallest boundary face containing z(j(sj)). This 
expresses that the normal component of the momentum may change, while the 
tangential component is conserved, when a light ray hits a boundary. 

Now one can describe the singularities of u using null bicharacteristics. Let o be 
the zero section of T*M. The location of the singularities is described by an object 

WF(«) C T*M\o = {(z,C) ■ Ct^O} 

that is conic in £, i.e. (z, £) 6 WF(tt) if and only if (z, rQ € WF(u) for every r > 0. 
WF(k) is called the wave front set of u, and it describes where (in z) and in which 
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codirection ( is the distribution u not C°° . More precisely, the definition of WF(u) 
is that (zo, Co) ^ WF(u) if and only if there exists (j) £ C£°(M), </>(zo) 7^ such that 
the Fourier transform J-{<pu) of <j)u is rapidly decreasing in an open cone around 
Co- Here we assume that M is boundaryless; otherwise we need to require that 
4> is supported in the interior of M. Again, there is a natural definition at dM 
which we do not give here. (There are more natural versions of this definition using 
pseudo-differential operators that I will describe later.) As an example, consider 
the step function: write z — (z%, z"), u(z) = 1 if z\ > 0, u(z) = if z% < 0. Then 

WF(«) - N*{ Zl = 0} \ o = {(0, z", Ci, 0) : Ci ^ 0}, 

the conormal bundle of the hypersurface z\ — 0, with its zero section removed. The 
same statement holds, with = possibly replaced by C, if we take any C°° function 
uq on M, and then define u — uq in z\ > and u = in z\ < 0. Informally, 
one might say that u is singular in z\ at z\ = 0, but it depends smoothly on z" , 
The wave front set thus pinpoints not only the locations z of singularities (lack of 
smoothness) in M, but it refines it by also giving the frequencies (or rather direction 
of frequencies) at which these appear at z. 

The theorem we are after is the following. In early versions it goes back to 
Lax |3j| , its boundaryless version is due to Hormander |28[], the smooth boundary 
versions are due to Melrose, Sjostrand, Taylor and Ivri i |33| , |43| , fil] , and the 
corner version in the analytic category is due to Lebeau J36[ (the C°° version is still 
not known in the corner setting) while a different extension, to conic points, is due 
to Melrose and Wunsch Ej|. 

Theorem 2.1. Suppose Pu £ C°°(M), and if dM ^ then u\ dM = 0. Then 
WF(u) C S = p~ 1 ({0}) (microlocal elliptic regularity). Moreover, WF(u) is a 
union of maximally extended generalized broken bicharacteristics inside £ (propa- 
gation of singularities). 

This theorem states that if a point (z, C) £ T*M \ o is in WF(u) and u solves 
Pu £ C°°(M), and satisfies a boundary condition if appropriate, then there is 
at least one maximally extended generalized broken bicharacteristic through (z, C) 
that is completely contained in WF(w). Of course, in the absence of boundaries, 
and often even in their presence, there is a unique maximally extended generalized 
broken bicharacteristic through (z, C), so the statement is that this bicharacteristic 
is completely in WF(«). However, as soon as codimension two or higher corners 
appear, there is no hope for such uniqueness, and this theorem is the optimal 
statement. 

At least in the nicest settings (no boundaries, or non-degeneracy assumption at 
the boundaries which are assumed to be smooth), this theorem can be improved 
significantly to predict not only the location, but also the amplitude of the singu- 
larities of u. 

3. Propagation in many-body scattering 

There is an analogous setup for scattering. Now we want to understand how 
interacting particles behave. Again, there is a classical mechanical setup (the ana- 
logue of geometric optics) and a quantum mechanical setup (the analogue of the 
wave equation). To focus on the most relevant points, I formulate the problem in 
a time-independent fashion, though it is easy to reformulate everything in a time 



dependent way. We only do this in a remark following Theorem 3.1 
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X23 

FIGURE 2. Collision planes X12, X13 and X23 and translates X' 12 
and X'12 ot X12. V\2 is constant along X12, is a (typically different) 
constant along X[ 2 , etc., so it does not decay at infinity unless it 
is identically zero. 

Thus, we want to understand tempered distributional solutions u of (H—X)u = 0; 
here A £ K. is the energy, and H is the Hamiltonian, i.e. the analogue of H — X 
is P above. Namely, if we have N particles, each of which is d-dimensional with 
positions x%, . . . , xn € R d , mass mi, . . . ,mjv, and the interaction between particle 
i and j is given by a potential Vij (which is a function on R d ), then the Hamiltonian 
describing this system is 

N 1 

H = 7^— A ^ + Vijfa - Xj) = A + V, 
t—i i<j 

which is an operator on (functions on) K™ = R^. Planck's constant h is here taken 
to be 1; it could be absorbed in the x\ by a simple rescaling. 

Now H is elliptic in the standard sense, namely its principal symbol is ^ 2m"l£l 2 ' 
which never vanishes outside the zero section o. Note that the potential is lower 
order than A in the standard sense, so it is not part of the principal symbol. So, 
by the previous theorem, 

(H - X)u = => WF(u) = N«e C°°(]R n ). 

So the only possibility of interesting behavior for u is at infinity, and this is exactly 
what we want to understand. 

The main feature of many-body problems is that even if Vij decays at infinity on 
M. d , it does not decay at infinity in K™ since it is a constant along Xij — {xi — Xj}, 
as well as along its translates X^ , X'/j , so it does not decay if we go to infinity, 
say, along X^; see Figure ||. The X^ are called collision planes (as are their 
intersections) since at X^ particles i and j are at the same place. 

In the two-body problem one actually has H — A X1 ^ X2 + Vi 2 (xi — x 2 ), i.e. V12 
still does not decay at infinity, e.g. if one keeps X\ — X2 but lets Xi — > 00. However, 
one can easily remove the center of mass by performing a Fourier transform along 
Xi2- This conjugates H — X to H 12 + |£i2| 2 ~ 'V where ^12 is the variable on 
Xl 2 , and H 12 = Axis + V12, X 12 being the orthocomplement of X12. Thus, one 
reduces the study of H ~ X to that of a Hamiltonian on X 12 , namely H 12 — A', 
A' = A — |£i2 1 2 being a shifted spectral parameter. Now V12 decays at infinity 
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Figure 3. On the left, broken geodesies in R n \{0}, n = 2, broken 
at the collision planes X a , Xb and X c . On the right, the projection 
of broken geodesies in R n \ {0}, n — 3, emanating from the north 
pole, to the unit sphere So, better understood as the sphere at 
infinity. The C , C& are the intersection of the collision planes X a , 
Xb with So; dimX a = 2, dim At, = 1. 



(we are working on X 12 !), so H 12 can be considered as a perturbation of Ajfi2, 
hence its analysis is rather simple. Notice that the point spectrum of H 12 gives rise 
to a branch of the continuous spectrum of H : this is a phenomenon that is very 
typical in many-body scattering. The center of mass can also be removed in any 
actual many-body problem, but one still obtains a Hamiltonian with non-decaying 
potentials as before. 

One can still talk about classical mechanics, just as for the wave equation, using 
bicharactcristics. These are deterministic - if V is smooth enough (we usually as- 
sume that V is C°°). But much like for corners, there is a compressed description 
of dynamics near infinity. This is somewhat more complicated than for the wave 
equation, but only because particles can be bound together. Thus, even the 'classi- 
cal' description is partly quantum. These two facts, the presence of collision planes 
and the bound states, are the two crucial features of many-body scattering. 

The compressed dynamics in the absence of bound states looks just like in the 
wave equation setting. One should think of this as a good description when a 
classical trajectory is uniformly near infinity. 

More precisely, it is convenient to introduce Agmon's generalization of the many- 
body problem, which amounts to using the vector space structure of R n as the 
setting. One can also give geometric generalizations (in the sense of differential 
geometry) that arose from the work of Melrose El], and I will do this later. 

So we work on the vector space Xo — R™, equipped with the Euclidean metric. 
We are also given a finite collection X — {X a : a e /} of linear subspaces X a of 
R", called the collision planes. We assume that X is closed under intersections, 
and X = R" e X, X t = {0} G X. We let X a = X^ be the orthocomplement 
of X a in R", so R™ = X a X a . (Agmon's generalization is thus that the X a do 
not have to come from intersections of the planes Xij = {xi — Xj}.) We write the 
corresponding coordinates as (x a , x a ), and denote the orthogonal projection to X a 
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by 7r a . A many-body Hamiltonian in potential scattering is an operator of the form 



where V a is a real valued function on X a in a certain class, for example V a is a 
symbol on X a of negative order: V a G S~ p (X a ), p > 0. We also assume that 
Vo = for normalization; note that X° = {0}, so Vo would simply play the role 
of the spectral parameter. We sometimes drop the pull-back notation from now on 
and write H = A + J2 a Va- 

Another useful piece of terminology is the following. We say that V a is short range 
if V a G S~ p (X a ) for some p > 1. We say that V a is long-range if V a G S~ p (X a ) for 
some p G (0, 1]. The Coulomb potential is thus 'marginally long-range', at least if 
we ignore its singularity at (which is not a serious problem anyway) . Whether V a 
is short- or long-range does not make any difference for the propagation phenomena 
we discuss in this section. However, it does make a major difference for the precise 
behavior of generalized eigcnfunctions at the 'radial sets' which we discuss later. 
This also shows up in the related issue of asymptotic completeness. 

Yet another notation we use on occasion is that of a fc-cluster. Physically, a 
cluster describes particles that are close (or collide), and a fc-cluster means that 
there are k clusters of particles, inside each of which the particles are close to each 
other. So in iV-particle scattering, the TV-cluster describes N asymptotically free 
particles (none is close to any other), hence we say that the collision plane Xq = W l 
is the iV-cluster. On the other hand, if X a ^ {0} is such that Xb C X a implies that 
X = {0}, then X a , or rather a, is a 2-cluster. E.g. given five particles, a 2-clustcr 
is where x\ = X2 and X3 — X4 = £5, i.e. the particles 1 and 2, resp. 3, 4 and 5, 
are close to each other. In general, a fc-cluster X a can be defined by the length of 
nested chains of collision planes inside X a . 

One need not assume that all interactions between the particles are via potentials. 
Indeed, V a may be allowed to be any first order differential operator on the vector 
space X a with symbolic coefficients of negative order. Also, one may generalize the 
metric g in an analogous fashion, as discussed later, which in effect allows V a to be 
second order provided that H remains elliptic. To simplify the notation, and due 
to the traditions, we mostly talk as if V a were potentials, but the generalization to 
such higher order perturbations requires only occasional and minor modifications, 
which will be pointed out. 

The subsystem Hamiltonians arc defined by 



Note that X a C X if and only if X a D X b , so above V& is really the pull-back of 
Vb from X b to to X a by the orthogonal projection. Thus, H a is an operator on 
(functions on) X a , and indeed it is a many-body Hamiltonian. 
We also let 



be the singular, respectively the regular, part of X a . Thus, if X c is a collision 
plane and X a is not a subset of X c , then X a n X c is a proper subset of X a , and 
is a collision plane (since X is closed under intersections), so X a C\ X c C A a:S i ng . 



H = A + ]T(7rTK 



a 



H a = A Xa + ]T V b . 



X a CXb 
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Correspondingly, V c decays at Ar a , r eg, so 

H a = A Xa +H a , 

which is an operator on (functions on) Xq = R™, has the property that H — H a 
is a function that decreases at X a ^ cg . So H a should be thought of as a good 
approximation of H at X ajTes . Note that X ajS i ng is a finite union of codimcnsion 
> 1 submanifolds of X ai so Af ajreg is in particular an open dense subset of X a . Also, 
note that Ax a plays a role analogous to the kinetic energy of the center of mass in 
the two-body setting, but now this description only valid locally, at X a>Tes . 

Having thus described the configuration space X = Xq = R", the next step is 
to describe the phase space, as was done first in f35| and The main goal in 

the process is to obtain a space on which broken bicharacteristics behave well. We 
remind the reader that we are concerned with singularities at infinity, hence with 
bicharacteristics that are uniformly close to infinity. Later we give a compactificd 
description, but here for simplicity we give its homogeneous version, much as for 
the wave equation where bicharacteristics were integral curves of the homogeneous 
principal symbol. So we start with T* X, but we wish to compress it at X a in such 
a way that at X a ^ cg , T^ a rcg X is replaced by T Xa re X a = T*X a ^ cg . For broken 
bicharacteristics this has the effect that only the X a -tangential component of the 
momentum is preserved at X aircg . So we define the compressed cotangent bundle 
as 

Note that this is at first just a set, equipped with a projection T* X — ► X \ {0} 
induced by the bundle projections T* X a rcg — > X ajTeg . There is also a natural 
R + -action on T*X via dilation in the configuration variables: 

(3.1) R+ x T*X a . Icg 3 (r,x a ,Q ~ (rx a ,^ a ) £ T*X a . Jcg . 

We topologize T*X via the projection 

n '■ T x\{o} X ~> T*X, 
whose restriction to X is the pull-back of one-forms by the inclusion map 

X a ,Teg X. Thus, writing (^ a , £ a ) as the momenta dual to (x a ,x a ), ir projects out 
the normal component of the momentum, £ a . The topology is then the weakest 
topology that makes it continuous, i.e. a set C in T'*X is closed if and only if ir^ 1 (C) 
is closed. 

We can now describe the contribution of the bound states to the characteristic 
sets. As mentioned above, this is one of the most interesting features of many-body 
scattering that has no analogue for the wave equation. These are conic subsets of 



T*X (conic with respect to the R+-action in (3.1)). The characteristic sets describe 
where certain operators are not elliptic, i.e. invertible, at infinity, in a precise sense 
described in the subsequent sections. They correspond to the 'energy shell', i.e. 
being on the characteristic set at energy A means that the particles have total 
energy A. We let 

Charo(A) = {(x,OeT*X: g{£) = A} 

be the free characteristic variety, with g being the metric function on T*X, and 
more generally we set 

Char a (A) = {(x a ,£ a ) e T*X a : \-g a {U) G spec pp (i? a )} C T*X a . 
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Char (A) n T p * X Char a (A) n T p * 2 X a Char 6 (A) n T p * 3 X b 

Figure 4. The characteristic set of many-body Hamiltonians. 
Here H is a 4-body Hamiltonian, a is a 3-cluster, b is a 2-cluster, 
Pi S Xo )re g 5 P2 S X aire g, P3 G ^b, r eg- The solid dots are the radial 
sets, denned below. 



Notice that A = g a (£,a) + f a , e Q G spec pp (7J°), corresponds to the splitting of the 
total energy A to the kinetic energy of the cluster, g a (^ a ), plus the energy of the 
bound state, e a . Thus, Char a (A) describes that particles may exist in a bound state 
of H a , of energy 6 Q , along X ai with kinetic energy g a 

(£a) = ^ — e a- Moreover, H° 
is the zero operator on X° = {0}, so if a = 0, these two definitions are consistent. 
If X a C Xb, the pull-back of one-forms gives a projection ir oa : T x X — > T*X a - 
Let 

Char(A) = U Char a (A) C f *X, 

Char a (A) = Ux b DA 7rf, a (Char b (A)) n T*X a , rcg . 

In order to understand Char(A) it is important to keep in mind several results 
on the structure of the eigenvalues of the subsystems. So let 

Aa = U b:X b £X a spec pp (H b ) 

be the set of thresholds of H a . Fundamental results of Perry, Sigal and Simon |33) 
and of Froese and Herbst jl6| show that A a is closed, countable, and the countable 
set spec pp (H a ) can only accummulate at A a , so 

A' a = A Q U spec pp (H a ) = U b:XbcXa spec pp (H b ) 

is also closed. Hence, Char(A) is a closed subset of T*X. In fact, the quotient 
E(A) of Char(A) by the R + action (which can be realized by restricting the various 
bundles to the unit sphere, Sq — {x € Xq : \x\ — 1}) is compact, and indeed it is 
metrizable, see Q. Since compact topological spaces have better properties than 
non-compact ones, it is quite natural to work with S(A), although we do not follow 
this route in this section. We also remark that it is much better to think of E(A) 
lying at the sphere at infinity, rather than at So, since it is the dynamics at infinity 
that is described here. We will take up this approach in later sections. 

We also recall another result of Froese and Herbst [jl5| , namely that eigenfunc- 
tions ip a of H a with eigenvalue e Q decay exponentially on X a , at a specified rate, if 
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e Q ^ A a . This generalizes to higher order perturbations, but requires a somewhat 
different approach, see |^|. In fact, this is the only place where second order per- 
turbations behave differently from first or zeroth order ones. For the latter, there 
can be no positive energy bound states, while for the former this has been only 
proved for small metric perturbations in p7[ , and it is not clear whether it holds 
more generally, especially for trapping perturbations. (Note that if H a , n / 1, is 
trapping then H is trapping at infinity!) 

A generalized broken bicharacteristic (at energy A) is then a continuous map 
7 : / — > Char(A), / an interval, such that a Hamilton vector field condition holds. 
To see what this is, we consider a subset of continuous functions ortT* X, namely the 
class of 7r-invariant C°° functions on T*X. 7r-invariance means that if £, C' <= T*X 
and 7r(C) = tt(C') then /(£) = /(£'). If / is 7r-invariant then it induces a function 
f n on T*X by f^{q) — /(C) if q = tt(C)- Moreover, if / is smooth (or indeed just 
continuous) then is continuous by the definition of the topology on T*X. 

Now, if 7 is a curve in a manifold, one way to put that it is an integral curve of 
a vector field V is that for all smooth functions /, 7)|s= So = {Vf)(^(so})- 
If / is a smooth it- invariant function on T*X, then / defines a C°° function on 
T*X a for all a, so H 9a f makes sense. Here H 9a is the Hamilton vector field of the 
metric function g a on T*X a , so explicitly, H 9a — 2£ a • d Xa . Now we would like to 
say that along a generalized broken bicharacteristic 7, 4-{f-K 7)| s =s should be 
given by H gb f for some b and some £ with 7r(£) = 7(s ). The problem is that there 
are many such points £ and clusters b, so this statement does not make any sense. 
However, we may replace the derivative by the liminf of the difference quotients, 
i.e. by D±h(so) = liminf s ^ So ^ s ]l^ s °^ , and demand an inequality instead of the 
equality. That is, we may demand that D±(f n o 7)(so) may not be less than the 
worst possible scenario as we run over all such b and £. Thus, the condition for a 
continuous map 7 : I — > Char(A) to be a generalized broken bicharacteristic is then 
that for any sq G /, if j(so) € T*X a TCg then 

D±(U o j)(s ) > M{(H g J)(Q : C e Char b (A), tt(0 = 7(^0), X a C X b }. 

If the set of bound states is discrete, then such a curve 7 is piecewise an integral 
curve of the Hamilton vector field of gb inside Char;, (A), where b may of course vary. 
In particular, if there are no bound states in any proper subsystem, the picture is 
very similar to wave propagation: the definition can be reduced to the analogue of 
Lebeau's f§. 

The structure of the generalized broken bicharacteristics, including the above 
claims, depends on having a large supply of 7r-invariant functions. But these exist, 
since the pull-backs of all functions on X to T* X is 7r- invariant, so one can localize 
in X using smooth cutoffs. Moreover, near X a ^ cg , each component of £ a is 7r- 
invariant, as is £ a • x a . Note that the generalized broken bicharacteristics depend 
on V, but only via the characteristic set Char(A), i.e. only via the bound states of 
the subsystem Hamiltonians. 

There is also a wave front set associated to many-body scattering which mea- 
sures the microlocal decay of tempered distributions at infinity. For a tempered 
distribution u, WF sc (u) is a closed conic subset of T*X. Apart from u, it depends 
on X , since T* X depends on X, but we suppress this in the notation, and write 

WF SC ( U ) = WF SC ,*( U ). 
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Its definition is slightly complicated, and I only refer to |6J| for the general defini- 
tion, which uses the structure of the pseudo-differential algebra, in particular the 
operator-valued nature of symbols at infinity. However, for generalized eigenfunc- 
tions of many-body Hamiltonians it is simple. Namely, suppose that (H — X)u S 
S(X), where S(X) is the space of Schwartz functions. For x = x a € X aties and 
£ a _€ X* we say that (x a ,£ a ) 4- WF sc (u) if there exists 4> G C^(X*) such that 
0(£a) 7^ and T~ x (f>Tu is rapidly decreasing in an open cone in X around x a . Two 
examples are: 

WF 8c (e^°) = 7r({(z,&) : x + 0}), £ E X*, 
WF^e™! 1 !) = n({(x, ^) : x ? 0}), a e R. 



•''I 

y„^i w> ^ i^oo^ 



More generally, if v is a symbol of any order on Xq, say u € S (Ao), and S C°°{Xq) 
is homogeneous degree 1 for |a;| > 1, then 

WF sc (e^v{x)) C 7r(graph#) = n({(x, (#)(z)) : M > !})• 
The condition |x| > 1 is due to the requirement of the homogeneity of <j> only for 
\x\ > 1; technically we should add a subset of |aj| < 1 to the right hand side to 
make it conic. The theorem on the propagation of singularities is the following. 

Theorem 3.1. Suppose that X € M and H is a many-body Hamiltonian. If 
u e 5'(M") and (H - \)u E S(R n ) then WF sc (u) C Char(A) (microfocal ellip- 
tic regularity). Moreover, WF sc (m) is a union of maximally extended generalized 
broken bicharacteristics of H — A (propagation of singularities). 

Remark. In the time dependent version, one considers tempered distributional so- 
lutions of (Dt + H)u = on Xq x R t . Then Dt + H still has the structure of a 
many-body Hamiltonian, with Dt + A in place of A, with collision planes given 
by X a x R, with {0} added for the sake of completeness. Thus, t is always a vari- 
able along the collision planes, so in particular, its dual variable r, is 7r-invariant. 
Moreover, Char a (A) is replaced by 

Char Q = {{x a ,t,i a ,r) e T*{X a x R) : -t - g a (£ a ) € spec pp {H a )} C T*X ai 

so effectively — r plays the role of the energy A. Generalized broken bicharacteristics 
can be defined as before with H gb replaced by H T+gb = d t + 2£ b • d Xb . The main 
additional issue is that they can only be expected to give a good description of 
propagation at finite ener gies since D t + H is not elliptic in the usual sense. So the 
analogue of Theorem [Tl] is that if it € S'(R" x R), (D t + H)u = 0, and u = ip(H)u 
for some ip £ C£°(R), then WF sc (w) is a subset of the characteristic set, and in fact 
WF sc (u) is a union of maximally extended generalized broken bicharacteristics of 
D t + H inside it. The proof of this statement only requires simple modifications of 
the proof of the theorem. 

The interpretation of the theorem is much analogous to that for the wave equa- 
tion. However, there is a difference which also occurs in the traditional microlocal 
setting for more general operators (i.e. for operators other than the wave operator), 
see |2^]. Namely, the orbits of the R + -action may be bicharacteristics, and then 
the statement of the theorem is empty at the points lying on these orbits since the 
wave front set is a priori conic. This happens for (x a , £ Q ) G T*X a)re g C T* X if and 
only if there exists some cluster b with Xb D X a , and £ £ Char^ (A) such that H 9b 
at C is tangent to the orbits of the R + -action. This happens, in turn, if and only if 
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£ a is parallel to x a and A — |£ a | 2 € spec pp (i? f> ). Such points are called radial points, 
and their collection is denoted by 

ft(A) = \J{(x a ,£ a )eT*X 

a,rcg ■ 3c G M, £ a CX a , 

a^l 

36, A b D X ai X — |£ a | 2 € spec pp {H b )}. 

As we discuss in Section ^, 7?.(A) plays an important role in asymptotic complete- 
ness. In many-body scattering it appeared in the work of Sigal and Soffer |55f| and 
was called 'propagation set' because in the time-dependent picture this is where 
particles end up as time goes to infinity. (In the stationary semiclassical picture, 
this is where non-trapped classical trajectories starting in a compact region end 
up.) It is thus unfortunate, in terms of terminology, that this is also the region 
where there is no real principal type propagation. 

Remark. In the time-dependent problem, the set of radial points is 

ft = lj{( x ^*'^' r ) G T *(^a,rcg x M) : x a = 2t£ a , 
a 

3b, X b D X a , -t - |£ a | 2 S spcc pp (H b )}. 

In terms of radial points, the difference between threshold energies A € A = Ao 
and non-threshold energies is that if A G A, then there are constant generalized 
broken bicharacteristics, i.e. bicharacteristics whose image is a single point. Namely, 
if x a £ A Qi i Cg and A £ spec pp (H a ), then 

(x a , 0) e Char Q (A) n T*A a , rog C Char(A), 

and H 9a = 2£ a ■ d Xa vanishes there, so (x a ,0) is indeed the image of a constant 
bicharacteristic. While this does not make any difference for the propagation of 
singularities, it does for the related limiting absorption principle, which in this 
generality is due to Perry, Sigal and Simon |p3| . 

Theorem 3.2. If X ^ A, then the limits R(X ± iO) = (H - (A ± iO))" 1 exists as 
bounded operators between L 2 a and H 2 _ s for s > ^. Here is the weighted Sobolev 
space (x)- l H m (R n ), l? s = H°. 

In fact, the proofs of the limiting absorption principle and the propagation of 
singularities are related. Indeed, the statement on propagation of singularities can 
be strengthened for R{X + iO)f 1 f G 5(E n ), by saying that WF sc (i?(A + iO)/) is not 
only a union of maximally extended generalized broken bicharacteristics, as follows 



from Theorem 3.1, but in fact it is a union of generalized broken bicharacteristics 
7 : R s -> f*X which go to 

ft+(A)=ft(A)n [J{(x a ,Z a ) e T*X a , ieg : x a •£>>()} 

a^l 

as s — > -co. That is, the singularities at 11+ (A) (where the statement of The- 



3.1 is empty) can only leave 1Z+ (A) in the forward direction. The limiting 



absorption principle is thus strengthened to: 

Theorem 3.3. If X (f. A, then for f E S(R n ), WF SC (R(X + i0)f) is a subset of the 
image of 1Z+{X) under the forward generalized broken bicharacteristic relation. A 
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similar statement holds for i?(A — iO)f with 1Z+(X) replaced by 

a) ^ T X a ^ re g . x a • £ a <c 0}, 

ajtl 

and the forward relation by the backward relation. 

In fact, if u G <S'(M n ) and WF sc (u) is disjoint from the image o/7£_(A) under 
the backward generalized broken bicharacteristic relation, then R(X + i0)u is defined 
by duality andWF sc (R(X + iO)u) is a subset of the image oflZ+(X)L)WF sc (u) under 
the forward relation. 

Remark. A ^ A can be also characterized by TZ(X) = 1Z+(X) UlZ- (A), i.e. that x a •£„ 
never vanishes on TZ(X) fl T*X a ^ og for any a. 

In the time-dependent setting, x a = 2££ a on 1Z, so x ■ £ a — It. So 1Z+, defined in 
1Z by x a ■ £ a > 0, is the subset of 1Z where t > 0. Hence the 'outgoing' terminology 
for R(X+iO) and 'incoming' for R(X— iO). In fact, the solution of (D t +H)u = with 
u| t=0 = 0, € S(Xo), say, is u(.,t) = e~ lHt <j). The time-dependent propagation 
of singularities shows that WF sc (u) is a subset of the union of the image of 1Z+ 
under the forward broken bicharacteristic relation and the image of 7Z— under 
the backward bicharacteristic relation. Using the spectral measure and Stone's 
theorem, 

u(.,t) = — [ e- lXt (R(X + iO)- R(X-iO)UdX 

Fixing some t/j E C^°(K), for <f> in the range of tp{H) we thus deduce that in t > 0, 
WF sc (u) arises from the R(X + iO) term, and in t < from the R(X — iO) term. 
So the time-dependent and stationary settings are very close: the only difference is 
that in the latter, A is a parameter, while in the former, it is a variable, A = — r. 

Again, one can make more precise propagation statements in some circumstances, 
such as three-body scattering, where the precise nature of the singularities can be 
analyzed, see |^2|, [|| . Here we only state the stronger implication for the structure 
of the scattering matrices, which we proceed to analyze. 



4. Scattering matrices 

Physically, the scattering matrices relate incoming and outgoing data in an ex- 
periment. In the time independent framework (where —A is the dual variable of 
time), for short-range potentials an incoming wave of energy A in channel a (a 
channel is the choice of a cluster a and an L 2 -eigenfunction ip a of H a of energy e a ) 
takes the following form in \x\ > 1: 

V = e- iy ^ l ^ l | 3;a |-^ i2 ^ i 3a,-( A)^(* a ) + W_ 

Fa| 

Similarly, an outgoing wave has the form 

u a ,+ = e i ^^^\x a \-"^^ 1 g a , + (^)^ a (x a ) +u' + 

pal 

i.e. the sign of the phase has changed. Here g a ,± may be taken e.g. L 2 functions 
on S a , the unit sphere in X a , or ideally, at least one of them may be taken C°° . 
In either case, u'± are 'lower order terms', namely they must be in L_ x , 2 . (Note 

that (x a )~ d '™ 1 * a 1 e L 2 s (X a ) for s < -1/2 but not for s = -1/2.) In fact, for 
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g a .± G C^ (S , ajreg ) we may take them to be of the form e~ l ^ x ~ ea > Xa > \x a \ 2 y 
where v is a Oth order symbol, with S aiieg denoting X ajlog n Sq. 

One can now produce tempered distribution with given incoming, or alternatively 
of given outgoing, asymptotics. A typical example is of the form 

(4.1) P a .+ {X)g a ^ = Uq ,_ - (H- (A + iO))-\(H - A)« a ,_); 

here the lower order terms may be dropped from without affecting u = 

Pa,+ (^)9a,— and g a ,- can be specified to be any smooth function on S a . In general, 



even if the incoming data are in a single channel a, as in (4.1), the corresponding 
generalized eigenfunction u of H will have outgoing waves in all channels. The 
S-matrix SqR(X) picks out the component in channel (3 by projection in a certain 
sense, see [p4[ . Thus, S a p{\) maps functions on S a , the unit sphere in X a , to 
functions on Sb, by 



for u as in (4.1). For example, the free-to- free (i.e. iV-cluster to TV-cluster in N- 
body scattering) S-matrix iSqo(A) maps functions on So, the unit sphere in K n , to 
functions on So, more precisely Soo(A) : L 2 (Sq) — > L 2 (So) is bounded. 

More precisely, let T+ be a pseudodifferential operator that is identically 1 on the 
outgoing radial set and identically on the incoming radial set; see the paragraph 



of (5.6) for a precise statement. Then 



(4.2) S a p(\) = i ((g-A)T + Pff,_(A)rP a , + (A), 

2i v /A - e/3 

i.e. for any g e C°°(5 a)reg ), ft G C°°(S b , Tes ), 

{h, S a p(\)g) = ((H - X)T + P^(X)h, i Pa,+(A)g). 

This is equivalent to the usual wave operator definition in the time-dependent 
setting, see An immediate consequence of the propagation of singularities and 
the definition of the scattering matrices is the following: 

Theorem 4.1. The wave front relation of S a p(\) is given by the broken bicharac- 
teristic relation. In particular, if no proper subsystem of H has bound states, the 
wave front relation of S'oo(A) is given by the broken geodesic flow on Sq at distance 

7T. 

While typically broken bicharacteristics can be continued in many ways when 
they hit a collision plane, it is important to keep in mind that under suitable 
assumptions (which rule out geometric complications) the broken bicharacteristic 
relation is Lagrangian, hence its dimension is the same as if there were no collision 
planes. The reason is that only a low dimensional family of broken bicharacteristics 
hits any specified collision plane, with the dimension of the possible continuations 
of each of these these bicharacteristics compensating to yield the correct dimension 
for Lagrangian submanifolds. 

The significant improvement in the three-body shown by Hassell and the 

author |33[ ^2j, is that one can pinpoint not only the location of the singularities, 
but also their precise form. This theorem was motivated by the geometric result 
of Melrose and Zworski |^6| , showing that the scattering matrix on asymptotically 
Euclidean manifolds is a Fourier integral operator. 
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Theorem 4.2. Suppose that H is a three-body Hamiltonian and the V a are Schwartz 
on X a for all a. Then iSoo(A) is a finite sum of Fourier integral operators (FIOs) 
associated to the broken geodesic relation on Sq to distance tt. Its canonical relation 
corresponds to the various collision patterns. The principal symbol of the term 
corresponding to a single collision at X a is given by, and in turn determines, the 
2-body S-matrix of H a at energies X' S (0, A). 

Remark. This result presumably extends to short range symbolic potentials, using 
the same methods, though it is technically more complicated to write down the 
argument in that case, and it has not been done. In fact, it should also extend 
to the A^-body problem, provided that there arc no bound states in any proper 
subsystem. Some assumption on the bound states is necessary, for otherwise the 
generalize broken bicharacteristic relation can become fairly complicated, see |36| . 
The reason why one does not need any assumption on bound states in three-body 
scattering is that for any 2-cluster a, Char a (A) D 7r 0o (Charo(A)) is either empty (if 
is not an eigenvalue of H a ) or consists of the boundary of 7To a (Charo(A)). In the 
former case there is no interaction (modulo smoothing terms) between the 0-cluster 
and the a-cluster dynamics, while in the latter case in the only place they interact, 
the two dynamics give the same propagation. 

It should also be noted that the normalization of S a p{\) is not the standard one 
in many-body scattering (which is based on wave operators) , but rather follows the 
geometric conventions |41j| . The difference is that in the wave operator approach, 
free motion is factored out, so the free scattering matrix is the identity opera- 
tor. On the other hand, in the geometric approach we describe the asymptotics 
of generalized eigenfunctions, or alternatively of the Schrodingcr equation. Since 
free particles move to infinity in the opposite direction from which they came, it is 
reasonable that the two should differ by (a constant multiple of) pull-back by the 
antipodal map, and this is indeed the case, see The distance 7r propagation 

along (not broken!) geodesies on the sphere indeed takes particles to the antipodal 
point. 

An immediate corollary, when combined with two-body results (e.g. analyticity 
of the S-matrix in A' and the Born approximation) is the following inverse result. 

Corollary. If the V a decay exponentially and dim X a > 2 for all a then Sqo(X) for 
a single value of A determines all interactions. 

This result is analogous to the recovery of cracks in a material by directing sound 
waves at it and observing the singularities of the reflected waves, except the last step 
which uses two-body results to get the potentials from the two-body S-matrices. 

The other extremal scattering matrices are the 2-cluster to 2-cluster ones, and 
they describe the physically most interesting events. Indeed, it is hard to make 
more than two particles collide in an accelerator, so the initial state in a physical 
experiment tends to be a 2-cluster. The following result is due to Skibsted p9| , 
and it also follows from the propagation of singularities and the definition of the 
S-matrices. 

Theorem 4.3. Let a and (3 be two-clusters, and suppose that either e a G spec d (_ff a ) 
and e/3 G spec d (H b ) , or V c is Schwartz for all c. Then the two-cluster to two-cluster 
S-matrix S a p(X) has C°° Schwartz kernel, except if a — (3 in which case the Scwartz 
kernel of S aa (X) is conormal to the graph of the antipodal map on S a , corresponding 
to free motion. 
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Thus, principal symbol calculations do not help in this inverse problem. Note 
that if H is a 3-body Hamiltonian, then e a 6 spec d (H a ) and E spec d (H b ) holds 
for any non-threshold bound state energies. The new result, in a joint project with 
Gunther Uhlmann, is the following JS2| . 

Theorem 4.4. Suppose that H is a 3-body Hamiltonian, a is a 2-cluster, a is a 
channel of energy e a <0, V a is a symbol of negative order (i.e. may be long range). 
For any p, > dim X a there exists S > such that the following holds. 

Suppose that sup \(l+\x b \) ti Vb(x b )\ < 6 for allb ^ a. Suppose also that I C (e a ,0) 
is a non-empty open set, and let 



R = 2 \J sup I — e a . 

Then S a > a "(X) given for all X E I and for all bound states a', a" of H a with 
e Q ',e Q " < sup/, determines the Fourier transform of the effective interaction V atC s 
in the ball of radius R centered at 0. 

The effective interaction is the interaction that arises if we consider the 3-body 
problem as a 2-body problem, i.e. if we regard the two particles forming the cluster 
a as a single particle. Mathematically, this amounts to projecting to the state 
ijj a in X a and obtaing a new Hamiltonian Aj o + V a , c s on X a . Thus, the effective 
interaction is physically relevant. Moreover, there is no hope for recovering anything 
better than V a<e s as shown by the high-energy inverse results of Enss and Weder 
0113), Novikov |9) and Wang 0, fnj. 



This theorem says that if the unknown interactions are small then the effective 
interaction can be determined from the knowledge of all S-matrices with incoming 
and outgoing data in the cluster a in the relevant energy range. In fact, near- 
forward information suffices as in two-body scattering, where this was observed 
recently by Novikov |5(| . Also, if one is willing to take small R and a is the ground 
state of H a , it suffices to know S aa (X) to recover V a ,eS in a small ball. 

In case V& decay exponentially on X b for all b ^ a, then V^eff decays exponen- 
tially on X a , hence its Fourier transform is analytic, so V a . c s itself can be recovered 
from these S-matrices. 

Remark. It is clear from the proof in [f32| that there is a natural extension of this 
theorem to many-body scattering at low energies. 

This result should extend to higher energies, i.e. sup / < is not expected to be 
essential. But it is hard to make R greater than 2^—e a even then. The reason is 
that our method relies on the construction of exponential solutions following Fad- 
deev Q , Calderon [|| , Sylvester and Uhlmann Q and Novikov and Khenkin Q , 
but in the three-body setting. One thus allows complex momenta p £ C(X a ), the 
complcxification of X a , and one wants to construct solutions of (H — X)u = of 
the form 

e^{i) a {x a )+v), 

where v = v p is supposed to be 'small' in the sense that it goes to as p — > oo 
in an appropriate fashion. Note that with v = these complex plane waves solve 
(H a - X)u = with 

(4.3) X = p-p + e a ; 

this expresses that the total energy A is the sum of the kinetic energy, p ■ p, and the 
potential energy e a . 
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To construct u, we need to find v, and its study reduces to that of the conjugated 
Hamiltonian 

e -ip.x a (jj _ x y P -x a =H a + A ^ + 2f) . D ^ +Ia ^ f a 



with p S C(X a ) the complex frequency. Here we used (4.3). Now, I a is considered 
as a perturbation (this is the reason for the smallness assumption in the theorem) , 
so we really study the model operator, 

H a + A Xa +2p-D Xa -e a . 

Taking the Fourier transform in the X a variables, one obtains 

H a + \i a \ 2 + 2p-^ a -e a . 

Writing p — zv + p± with \u\ = 1, px • v = 0, p, v real, zeC, this operator becomes 

H a + |£ a | 2 + 2 P± ■ £ a + 2zu -L-t a = H a + (£, + Pl _) 2 + 2zv ■ £ a - \ Pl _\ 2 - e a . 

If p is not real, then neither is z, so this operator is invertible if v ■ £ a ^ since H a 
is self-adjoint. On the other hand, if v ■ £ a = 0, this operator becomes 

H a + (da + p ± j 2 -\p ± \ 2 ~e a , 

i.e. its invertibility properties correspond to the behavior of the boundary values 
of the resolvent of H a at the real axis. If |pj_| 2 + e a < 0, i.e. if \p±\ < y/—e^, 
then the spectral parameter \p±\ + e Q — 

(Co + P±) 2 is negative, so only the bound 
states of H a contribute to the characteristic variety, i.e. the two-cluster a may not 
break up. On the other hand, if \p\ > ^/—t ai such a break-up is possible even if 
A < 0, i.e. where the break up may not happen for real frequencies. The break-up 
greatly influences analyticity properties, hence one cannot easily use large p±. On 
the other hand, one needs such large pj_ to recover V a , e f[ on larger balls, hence the 
limitation in the theorem. This also suggests that the fixed energy problem would 
be hard, since then one always needs to let p± — > oo to keep p- p = \p±\ 2 + z 2 fixed 
and yet have p — > oo. 

5. Many-body scattering pseudo-differential operators 

I will present the calculus from the compactified point of view. Both the one- 
step polyhomogeneous (i.e. 'classical') and the non-polyhomogeneous calculus can 
be described in non-compact terms, i.e. directly on Xq, but this is more complicated 
and less natural. Indeed, one of the beauties of compactification is that it exactly 
captures the structure of many-body Hamiltonians. We warn the reader here that 
from now on the Euclidean variable is written as z, rather than x in the preceeding 
sections, for compatibility with previous papers espousing this approach, such as 



To see how the compactification should go, recall first that a classical symbol of 
order on R™ has an asymptotic expansion 

oo 

a{ru) ^Y J ^ j a j (oj), dj E C^S™" 1 ), 
j=o 

in the polar coordinates (r, to): z — ruj. The meaning of such an expansion is that, 
for any k, the difference of a and the sum of the first k terms on the right hand side 
is a symbol of order —k. This expansion is just a Taylor series at r — oo, or rather 
at V -1 = 0'. So we compactify R™ into a ball B n by adding points (0, w), u) € §" _1 , 
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and making (r -1 ,!^) = {x,uj) coordinates near these points. The resulting space is 
called the radial compactification E™ of R". Thus, a classical symbol of order is 
simply a smooth function of 1"; the asymptotic expansion at infinity is its Taylor 
series around the boundary, x = 0. 

This compactification, whose utility in this context was emphasized by Melrose 
|4l|| , can also be realized as the closed unit upper hemisphere via a modified stere- 
ographic projection. So let 

RC : R" - §1, RC(z) = {-^, -^-), (z) = (1 + \z\ 2 f' 2 , z e R". 

Then n of the n + 1 variables (tjtj 7§y) give local coordinates on various regions of 
In particular, in coordinate patches near the equator, which is <9§" , (or 
indeed x = and n — 1 of (or indeed uij = j^) can be taken as coordinates, 

showing that §" can be identified with the radial compactification R n . A slightly 
modified version of x (it needs to be smoothed at z = 0, where 'x — oo'), or 
can be taken as a boundary defining function. We will usually write x for this, 
so x = I I 1 for \z\ > 1, say. (A boundary defining function is a non- negative 
function whose zero set is exactly the boundary, and whose differential does not 
vanish there.) 

How can we adapt this to many-body scattering? Let X a denote the closure of 
X a in the compactification R™ of R™, and let C a = dX a C <98™ = Co- The closure 
of any translate of X a intersects Co in the same submanifold (a sphere) as X a itself. 
Indeed, writing the coordinates as (z ai z a ) on Xq = X a ® X a , local coordinates near 
C a are given by Z a = 1 and dimA a - 1 of Thus, Z a -> as x -> 

along any translate, since z a is constant along these. So V a is not even continuous 
on Ao, as it takes different values on the different translates of X a . However, it is a 
negative order symbol (in particular continuous with boundary value 0) on Ao\C a , 
if V a is such onl"; see Figure |[ 

So the compactification works for V a , except at C a - To remedy this, we blow 
up C a - This is an invariant way of introducing polar coordinates about it (i.e. 
projective coordinates in various charts). That is, curves approaching C a from 
various normal directions will correspond to different points on the blown-up space 
[Xo;C a ]- Since C a is given by x = 0, Z a = 0, in local coordinates, this means 
concretely that the components of Z a jx become coordinate functions on the part 
of [X ; C a ] where this quotient is finite. (For the sake of completeness, a complete 
set of coordinates in this region is given by x, the components of Z a as well as the 
dim X a — 1 coordinates on the sphere y a = rf^r! see Figure ^.) But Z a /x = z a , so 
it is now easy to see that for classical symbols V a on X a (of negative integer order), 
V a is a C°° function on [A ; C a ]. 

In general, there are many collision planes, and we blow them up recursively, 
starting with ones of the largest codimension, to get [Aq;C], 

C = {C a : X a e X, a ? 1}. 

We refer to Q for details. 

There is no reason at all to take Ao as the space we start with. Given any 
compact manifold with boundary, A, and a cleanly intersecting family of closed 
embedded submanifolds C of dX, we can define [X;C] analogously. For instance, 
one can start with A = R™ x S fc . The space [A;C] is equipped with boundary 
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Figure 5. Translates of X a on [Xo;C a ]. 




Figure 6. The blow up of C a , given by Z a = 0, x = 0. 



fibrations given by the blow-down maps, see |58|] for a simpler case where these first 
appeared explicitly. 

Having described the configuration space, we turn to differential operators. X$ 
has a nice algebra of differential operators, consisting of operators with symbolic 
coefficients: J2\ a \<m a a(z)D^ , a a £ S°(Xq). We may require instead that a a is 
'classical', i.e. that a a G C°°(Xq). The resulting algebras were denoted Diff scc (Xo) 
and Diff sc (Xo) by Melrose; he called them 'scattering differential operators'. 

This setup generalizes to the geometric set-up as follows. Let (x,y), y = 
(yi, . . . ,y n -i), be local coordinates near dX . Then the vector fields in Diff sc (X ) 
are linear combinations of x 2 D x and the xD yj with coefficients in C°°(Xq), as can 
be seen easily by an explicit calculation. 

Now, if X is a manifold with boundary, Vb(X) is the Lie algebra of vector 
fields tangent to dX, and V sc (^) = xVb(X), where x is a defining function of dX. 
V sc (^0 is independent of the choice of x. Then Vb(^0 is spanned by xd x and d y 
over C°°(X), so V sc (^) is spanned by x 2 d x and xd y over C°°(X). By definition, 
these generate Diff sc (X). Also, V sc (^) is the set of all smooth sections of a vector 
bundle over X, this is denoted by SC TX. Its dual bundle is the scattering cotangent 
bundle, denoted by SC T*X. In the Euclidean setting, 

sc TXo = Xo x Xq, sc T Xq = Xq x Xq. 
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The way to generalize this differential operator algebra to one that includes 
many-body potentials is to allow singular coefficients a a £ C°°([X;C]). Thus, 

Diff sc (X;C) =C co ([X;C])® c ^ { x ) Diff sc (X). 

In particular, if if is a many-body Hamiltonian, with either potential or higher 
order interactions, then H £ Diffg C (Ao; C). 

We let SC T[X;C] = f3* sc TX and SC T*[X;C] = (3* SC T*X, where f3 : [X;C] -> X 
is the blow-down map, and we are pulling back the vector bundles by it. Hence in 
the Euclidean setting, 

SC T[A ;C] = [X ;C] x X , SC T*[X ;C] = [X ;C] x X*. 

Now it is natural to define pseudo-differential operators using these bundles. Al- 
though I restrict the discussion to the Euclidean setting, the construction general- 
izes to any X via localization. 
So we consider symbols 

(5.1) ae(z)- l (C) m C°°([Xo;C]xX*), 

Xq = R™, where ( is the dual variable of z, i.e. the variable on Xq. Note that this 
means that a is a classical symbol of order m in £. As usual, we dehne the Schwartz 
kernel of the left quantization of a by 

(5.2) A = q L (a) = (2vr)-" J e^-^<a(z, dC, 
understood as an oscillatory integral. In particular, for any / £ <S(M n ), 

Af(z) = (2tt)- b / / e^-^<a(z,C)f(z')dCdz', 



again understood as an oscillatory integral. We write A £ fy sc (Xo;C) for this class 
of operators. 

Note that (z) l a £ S™(Xo]Xq), Hormander's uniform symbol space |^7j Sec- 
tion 18.1], so A — (z) l A, A £ <f™(Ao), the uniform ps.d.o. -algebra arising by 
quantizing S™(Xq]Xq) as in ([5.2]). In particular, since the mapping properties of 
^^(Ao) between weighted Sobolev spaces H r ' s are well known, the corresponding 
properties of A follow. Namely, A : H r ' s — > jj r — m > s + l for all r, s, where 

H r ' s = (z)- s H r = (f £ S'QfT) : (z) s f £ H r }. 

Now, ^ sc (Xq;C) is a *-algebra, in particular is closed under composition. Indeed, 
since x l/ sc (Ao;C) C x P 00 (Ao), and the latter is closed under composition, it suffices 
to follow the usual proof and make sure that the product is in ^ sc (Xq;C), rather 
than merely in ^^(Xo). Thus, the key fact is that for any 

b£ {z)- l {Q m C°°([Xo;C\, x [X ;C] Z , x (A*) c ) 



there exists a as in (5.1) such that the induced operators 



(5.3) B = {2n)- n e l{z ~ z "><b(z 7 z',()d(, 



and A as in (5.2) are the same. The proof of this claim is standard. Indeed, we 
can expand b in Taylor series in z' around z = z' to finite order k. The finite order 
terms depend on z 1 only via [z! — z) a , \a\ < k. We rewrite (z' — z) a e l ^ z ~ z ''^ as 
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(— l)l Q l£)"e^ z z and integrate by parts with respect to £. Thus, the a-term is 
the left quantization of 

(5.4) -D^,D^b(z,z',OU'=z, 

which is of the desired form, i.e. is in {z)~ l {Q) m C°° {[Xq]C] xXq). In fact, the weight 
(C) m can be replaced by (£) m ~H due to the symbolic properties of b in (Xg)f, but 
no corresponding change may be made for the z weight. Similarly, the remainder 
term is of the form 

K k {z,z') = {2-K)- n [ e^-^<b k {z,z',Qd(, 

(5.5) JK" 

b k G (z)-'(C) m -' £ - 1 C 00 ([Xo;C] 2 x [X ;C] Z , x (X*) c ). 
Now we can asymptotically sum the b a to get a new symbol 

ce(.)- ! {crr([x ;C]xX *). 



Let C be the left quantization of c. Then B — C is of the form (5.5) for all k, with 
bk replaced by some b' k with the same properties. It is then straightforward to show 
that the Schwartz kernel K' of B — C is C°°, decays rapidly with all derivatives as 
(z — z') — > oo, and more precisely it is of the form 

if' eC°°{[X ;C] z x (* )*-«0 

with infinite order vanishing at the boundary of the second factor. Taking its 
Fourier transform b' in z — z' , K' is thus the left quantization of a = c + b', proving 
the claim, hence in turn that ^ SC (X ;C) is closed under composition. 

In the two-body setting, where C = 0, there is a principal symbol at infinity. 
Namely, if A G ^™' l (X), A = (^(a), then cr m ,i(A) is given by the restriction of 



[z 



(C)-' m a G C°°(X Q x X*) to d(X a x X*) = (dX Q x_X*) U (X x 9X*). Of the 
two boundary hypersurfaces, the restriction to Xq x dX^ yields the usual principal 
symbol, while the restriction to OXq x X^ is the principal symbol at infinity. More 
precisely, if I = 0, we can indeed define the part of a m fi(A) at infinity to be 
the restriction of a to (OXq) x Xq. The principal symbol is multiplicative, i.e. 
a m+m ,, lw (AB) = a m ,i(A)a m >,i f (B). Thus, [A,B] G * s ™ +m '- 1 ' z+z ' +1 (X), and its 
principal symbol is given by the Poisson bracket of their symbols, see Section |^. 



Since in the many-body setting we do not gain decay in z in (5.4), we cannot 
expect to have a commutative principal symbol at infinity, i.e. at 8[Xq]C]. For 
C G *2'°(X,C), y a G C Q , rcg , (a G X*, we let 

6„(y„,Ca) = (2^)- dimX ° / e^'-^'^'cd/a.z-.OdC G S'(X a x X a ) 



be the operator valued principal symbol of C at (y a ,Ca)- Thus, C a {y a ,C >a ) is a 
tempered distribution on X a x X a (denoted by the variables (z a , (z a )')), and it is 
in fact a many-body ps.d.o. itself: C a (y a ,C) G ^'"'''(X , C a ) corresponding to the 
collision planes X a n Xj, with & satisfying X, 3 X a . We also call it the indicial 
operator of C to make it clear we are not talking about the standard principal 
symbol. We also write C a (z a ,(a) in the same setting, where we extend C a (y a ,(a) 
to be homogeneous degree in z a . It can be easily seen to satisfy 

A a B a = (AB) a , 
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where on the left hand side we compose the operators A a (z a ,^ a ) and B a (z a , Ca)- 
Thus, multiplication of operators is only partially commutative, even to top order. 
This can be observed already from [D Za , V a ] = 0, hence certainly lower order at 
infinity, while [D z a,V a ] £ C°° ([Xq;C]) without any decay at C a - 

This observation has important implications for the positive commutator esti- 
mates that we take up in the next section. Namely, H must commute to leading 
order with the operators we want to microlocalize with. This means that these oper- 
ators A must have A a commute with H a , and the most reasonable way of achieving 
this is to have A a be a scalar multiple of ip(H a ), where e.g. tp E C^°(R). This mul- 
tiple defines a function on T*Xq] we want this to arise from a smooth 7r-invariant 
function for our estimates. On the other hand, ip(H a ) provides localization at the 
characteristic set. 

Here, however, I would like to talk about pseudo-differential constructions first. 
Namely, if A ^ R, or indeed if A G C \ [inf A, +00) then there exists a parametrix 
G(A) e *- 2 '°(X ;C) for H - A, i.e. such that 

(H - A)G(A) - Id, G(X)(H - A) - Id e ^-'"''"(XojC). 

Then the parametrix identities show that 

A S C\spec(ff) {H-X)- 1 S ®- 2 <°(X ;C). 

The parametrix construction proceeds inductively by constructing (H a — A) -1 in 
ty~ c ' (X a ,C a ) for every a / 1 and then combining these: there exists a Go(A) e 
ty~ c ' (Xq;C) with specified indicial operators (H a — A) -1 , hence satisfying 

(H - A)Gq(A) - Id, G (X)(H - A) - Id e f" 1 ' 1 ^^). 

Then the standard Neumann series argument yields G(A). 

The Helffer-Sjostrand argument |24| then shows that for any <fi £ C£°(R), 

<t>(H) = / d x 4>(X)(H - A)" 1 dX A dX, 

where <p is an almost analytic extension of 0: <\> € C£°(C), \d\4\ < Ck \ Im A| fe for all 
k. We can control (H — A) -1 in ' (X , C) as A — > R with semi-norm estimates 
0(| Im Al"- 7 ) (j depends on the norm), so we conclude that cf>(H) e ^°°'°(X;C). 

We can now explain the precise specifications on T + in ( |4.2| ) . Namely, we require 
that on a neighborhood of 7?. + (A) in T*X, the indicial operators T + should equal 

<fi(H) for some <fi £ C^°(R) identically 1 near A, and on a neighborhood of 1Z-(X) 
they should vanish. Explicitly this can be arranged by taking any </> as above, and 
any x £ C°°(R) identically 1 on (^,+00), identically on (-00,-^). Then 
let T + — (j)(H)qn(x{j§j)) , with qn denoting the 'right quantization' (i.e. where we 
take b = x(fpy) in (|^))- Although q R (x(^t)) is not in ty sc (X;C), due to the 
non-symbolic behavior of b as C - * °°, T+ is, namely T + 6 ^~°°'°(A > ; C), since 
</>(i?) is smoothing: see pE[ . Moreover, 

(5-6) 7^ a (z a , Ca) - X ( f fr i )^a(^, Ca)), 

hence has the desired properties. 

Our construction of 4>(H) in fact shows that if all potentials are in S~ p (X a ), 
p > 0, and Xa £ C°°(Xo) is supported away from G such that G& D C a does raoi 
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hold, then % a (0(H) - <j){H a )) £ ^ SC °°' P (X; C), hence trace class if p > dimX . In 
the three-body setting this shows that 

0(H) - 0(H O ) - ]T ^ H b) ~ ^ H *)) 

#b=2 

is trace class. Indeed, near C a this can be written as 

(0(H) - 0(H O )) - E (^)-W, 

and now all terms in parantheses are in 1 $'^ ! DO ' / '(X;C) near C a . So we conclude, 
with a proof that shows much more, a result of Buslaev and Merkureev: 

a{4>) = tr((0(H) - <f>(H )) - £ (0(H 6 ) - <f>(H ))) 

#6=2 

defines a distribution a £ C _00 (R). Writing a — defines the spectral shift 
function, up to a constant, which in turn, in two-body scattering, is the well- 
known generalization of the eigenvalue counting function on compact manifolds. 
These statements, as well as the following theorem, which is joint work with Xue- 
Ping Wang J68| , generalize to arbitrary many-body Hamiltonians (with short-range 
interactions as indicated). 

Theorem 5.1. Suppose H is a three-body Hamiltonian with Schwartz interactions: 
V a £ S(X a ), and that all interactions are pair interactions (i.e. V a 7^ implies that 
a is a 2-cluster.) Then the spectral shift function a is C°° in R \ (A U spec pp (H)), 
and it is a classical symbol at infinity (i.e. outside a compact set) with a complete 
asymptotic expansion: 




Note that a decays one order faster than in 2-body scattering, and two or- 
ders faster than Weyl's law on compact manifolds. This is because (j>(Ho) + 
S#&=2(^(H6) — 4>(Hq)) is, in a high-energy sense, closer to 0(H) than <j)(Ho) is 
to 0(H) in two body scattering. If not all interactions are pair interactions, the 
order of the leading term changes, namely becomes A^~ 2 as in 2-body scattering. 

The proof of this theorem relies on the propagation of singularities, applied to the 
Schwartz kernel of the resolvent, R(X + iO). (In fact, the theorem should generalize 
to symbolic potentials, but the proof would require a more precise microlocalization 
than provided by WF SC .) So we now turn to the positive commutator methods that 
prove this. 

6. MlCROLOCAL POSITIVE COMMUTATOR ESTIMATES 

First I sketch, somewhat vaguely, the idea of positive commutator estimates. So 
suppose that we want to obtain estimates on the solutions of Pu = f, where / 
is known, and is 'nice', and P is self-adjoint. Suppose that we can construct an 
operator A which is self-adjoint and is such that 



(6.1) 



i[A,P] = B*B + E. 
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Here B*B is the positive term, giving the name to the estimate. The point is that 
we can estimate Bu in terms of Eu. Indeed, at least formally, 

(u,i[A,P]u) = (u,B*Bu) + (u,Eu) = \\Bu\\ 2 + (u,Eu). 

On the other hand, 

(u,i[A,P]u) = {u,iAPu) - (u,iPAu) = (u,iAPu) + (iAPu,u) = 2Re(u,iAPu). 
Combining these yields 

(6.2) ||Bu|| 2 < 2\Re(u,iAPu)\ + \{u,Eu)\. 

This means that Bu can be estimated in terms of Pu, which is known from the 
PDE, and Eu, on which we need to make assumptions. The typical application is 
that E is supported in one region of phase space and B in another, in which case 
we can propagate estimates of u. 

In fact, one can also apply this estimate if one does not know a priori that 
Bu G L 2 . Namely, an approximation argument gives that if Pu and Eu are in 
appropriate spaces so that the right hand side of (|6 . 2h makes sense, then Bu S L 2 , 



and (3.5) holds. Considering pseudo-differential operators A of various orders, this 
means that we obtain microlocal weighted Sobolev estimates for u. Also, typically 
one has an error term F, i.e. i[A, P] = B*B + E + F, but F is 'low er order' in 
some sense. Thus, \(u,Fu)\ is added to the right hand side of (|6.2| ), but being 
'lower order' means that \(u,Fu)\ automatically makes sense, hence is irrelevant 
when proving that Bu € L 2 . 

In fact, this method also yields estimates for the resolvent very directly. Since 
for t e K, i[A, P-it] = i[A, P], and 

(u, i[A, P - it]u) = (u, iA(P - it)u) - (u, i(P - it)Au) 

= (u,iA(P-it)u) + (iA(P + it)u,u) = 2Rc(u,iA(P - it)u) - 2t(Au,u). 

Thus, we deduce 

(6.3) \\Bu\\ 2 + 2t(Au,u) < 2\Re(u,iAPu) \ + \(u,Eu)\. 

This is in particular an estimate for \\Bu\\ provided that t > and A is positive. 
Here we may take u — ut = (P — it)~ l f , defined for t > 0, say, and we find a 
uniform estimate for But as t — > 0. 

The question is thus how one can produce operators A which have a positive 
commutator with P as above. First, we recall how this happens in the scattering 
calculus. Namely, if A G ^T/(X), P e (X) then [A,P] € y™+ m '~ 1 , l + l '+\x) 

and 

0Vn+ m ' +i(i[A,P]) = H a p = -H p a, a = a m .i(A), p = <7 m ,. v (P), 

where H a is the Hamilton vector field of a, H p the Hamilton vector field of p. So 
modulo lower terms, which I ignore here and which are easy to deal with, we need 
to arrange that 

(6.4) H p a = -b 2 + e, 

and then take B, E with u(B) — b, u(E) = e. Indeed, under these assumptions ( |6.2[ ) 
shows that \\Bu\\ can be estimated in terms of Pu and Eu. That is, u microlocally on 
supp b is estimated by u on supp e (and Pu) in the precise sense described in the next 
paragraph, so we can propagate estimates of u from supp e to supp b. (Incidentally, 
this is a good example of the F term: only the principal symbols of E and B are 
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specified. Take any E and B with these principal symbols, F — i[A,P] — B*B — E 
is lower order.) 

This can be used in a very straightforward manner to obtain bounds on WF sc (w). 
Namely, one works with 'relative wave front sets', relative to x s H r = H r ' s , that 
is. Thus, for X — R", (z, £) ^ WFgj. s (u) means that there is a cutoff function <j> G 
C£°(R") with 4>(() ^ such that J 7 ^ 1 cj>!Fu is in H r,s in an open cone around z. But 
this is equivalent to the existence of some Q G W'^r'tX) such that a(Q)(z, C) 7^ 
and Qu G L 2 . Note that a(Q)(z, £) 7^ means that Q is elliptic at [z, Q. So if 

we find A G and consequently B G ^^^{X), E £ *™~ U+1 (X) as 

above, then the conclusion is that (if Pu G C°°(X)) 

WF^T'^ (u) nsuppe = =S> WFs? 1,_1 * i (tt) n suppb = 0. 

In scattering theory m is usually irrelevant by standard elliptic regularity. Thus, 
one iteratively reduces I, proving that supp b is disjoint from the wave front set with 
respect to more and more decaying Sobolev spaces. (In fact, b is shrunk slightly 
during the iterative procedure for technical reasons.) 

I now illustrate how to prove the propagation of singularities at OX for real 
principal type P € ^^°(X). For example, we may take P = H — A, A > 0, 
and H is a two-body Hamiltonian. (Note that microlocal elliptic regularity is the 
consequence of the standard microlocal parametrix construction.) We thus want 
to prove that if Pu G C°°(X) (or a microlocal version of it holds), f G SC T*^X 
and there is a point on the backward bicharacteristic through £ which is not in 
WF sc (u), then £ ^ WF sc (u). In fact, by a simple argument it suffices to prove that 
there exists a neighborhood U of £ such that if there is a point £ in U which is 
also on the backward bicharacteristic through £ and which is not in WF SC («), then 

eg- wf sc (u). 

The standard proof proceeds via linearization of H p , see pg|| . Thus, first note 
that x~ 1 H p is a smooth vector field on SC T*X which is tangent to the boundary. 
(For example, for Euclidean two-body Hamiltonians, 

x~ 1 H p = 2\z\(-d z = 2Y,Qz\d Z] , 

3 

and \z\d Zj is a smooth vector field tangent to dX, i.e. it is in Vh(X).) Thus, given 
any point £ G SC T*^X one can introduce local coordinates (q%, . . . , (fen-l) = feb q") 
on SC T*^.X centered at £ such that x~ l H p — d qi at SC T*^.X. Thus, bicharacteristics 
at dX are curves q" = constant. Now let X i S C™(R qi ) and X 2 € C^°(R^™ _2 ) be 
smooth functions supported near with the property that 

Xi = -bl + ei, 

61, ei G C£°(R), and suppei C (—00, 0). Let 

a = X1X2, b = 61X2, e = e x x\- 



Then (6^) holds. In fact, we can even allow weights and take 



x s Xixl,s£ 



since (x 1 H p x s )xi can be absorbed in x"(x 1 H P xi) by choosing x'x large compared 
to xi- This gives microlocal weighted L 2 estimates in a; _s_1 / 2 i 2 . The iterative 



26 



ANDRAS VASY 



i 


q" 




































ppa 































Figure 7. suppo superimposed on the linearized Hamilton flow, 
suppe is the shaded region on the left. 



argument, in which we gradually let s — > — oo, then allows one to conclude that 

suppe n WF sc (m) = and suppa n WF sc (Pu) = {b > 0} n WF sc (u) = 0. 

By choosing suppxi and suppx2 appropriately, we may arrange that e is supported 
near £ so that suppe n WF sc (u) = 0, and so that &(£) > 0, as shown below. 

There are several directions from here. One can use finer notion of regularity, 
such as Lagrangian regularity, which would correspond to using X2 that vanishes 
simply on a Lagrangian submanifold, or such as regularity at radial points, which 
is the subject of a joint paper with Andrew Hassell and Richard Melrose [ p3[ . 

Here I will talk about a rougher version, namely what happens if the bicharac- 
teristic 'flow' is more complicated, e.g. in the presence of boundaries or corners for 
the wave equation [ [43| , |36| or many-body scattering. In fact, here I will not explain 
the detailed behavior of bicharacteristics, rather just show how to microlocalizc 
positive commutator estimates in a versatile fashion. This method goes back to the 
work of Melrose and Sjostrand Q. 

The main point is that if we cannot put the operator P, or at least its Hamilton 
vector field H p in a model form, the previous construction will not work. Indeed, 
unless H p X2 = 0, H p (xiX2) will always yield a term XiH p X2y which cannot be 
controlled by {H p xi)X2'- the problem being near the boundary of supp%2- So 
instead use a different form of localization. First let r\ 6 C°° ( SC T* X) be a function 
with 

= 0, H p r,(0 > 0. 

Thus, rj measures propagation along bicharacteristics, e.g. r\ = q\ in the above 
example would work, but so would many other choices. We will use a function ui 
to localize near putative bicharacteristics. This statement is deliberately vague; at 
first we only assume that u> € C°° ( SC T* X) is the sum of the squares of C°° functions 
<jj, j — 1, . . . , I, with non-zero differentials at £ such that dq and daj, j = 1, . . . , I, 
span T^ sc T*j^X . Such a function to is non-negative and it vanishes quadratically 
at £, i.e. = and du>(£) = 0. An example is lo = q\ + . . . + q\ n _\ with the 
notation from before, but again there are many other possible choices. We now 
consider a family symbols, parameterized by constants 5 6 (0,1), e € (0,<5], of the 
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Figure 8. suppa in (j), a) coordinates, suppe is again the shaded 
region on the left. 

form 

a = Xo(2 Xi( — j- + 1 , 

e eo 

where 

(j> = rj + -to, 
e 

Xo(t) = if t < 0, X o{t) = e" 1 /* if t > 0, xi G C°°(R), su PPXl C [0,+oc), 
suppx'j C [0, 1]. Although we do not do it explicitly here, weights such as x s can 
be accommodated for any s £ R, by replacing the factor Xo(2— 7) by y (^4 ( 7 1 (2— 7)) 
and taking Aq > large. 

We analyze the properties of a step by step. First, note that </>(£) = 0, H p (f)(£) = 
HpViO > 0, and XiC^-jr + 1) is identically 1 near £, so H p a(£) < 0. Thus, H p a has 
the correct sign, and is in particular non-zero, at £. 

Next, 

£ S suppa =4> 0(£) < 2e and 7/(£) > -<5 - eS. 

Since e < 1, we deduce that in fact 77 = > — 25. But a; > 0, so <\> — 0(£) < 2e 
implies that 77 = — e _1 oj < <fi < 2e. Hence, w = = e(</> — 77) < 4e<5. Since a; 
vanishes quadratically at £, it is useful to rewrite the estimate as lo 1 / 2 < 2{e8) 1 / 2 . 
Combining these, we have seen that on supp a, 

(6.5) -5 - eS < 77 < 2e and uj 1/2 < 2(e5) 1/2 . 

Moreover, on supp a n supp x[, 

-S-eS < 77 < -<5 and cj 1 / 2 < 2(e<5) 1 / 2 . 

Note that given any neighborhood U of £, we can thus make a supported in U by 
choosing e and 5 sufficiently small. Below we illustrate the parabola shaped region 
given by supp a. 

Note that as e — > 0, but 6 fixed, the parabola becomes very sharply localized at 
lo = 0. In particular, for very small e > we obtain a picture quite analogous to 
letting suppx2 — ► {0} in Figure [?]. 
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So we have shown that a is supported near £. Wc define 

e = X o(2-^)fipXi(^ + l), 

so the crucial question is whether -ff p </> > on supp a. Note that choosing So € 
(0, 1) and eo G (0, <5o) sufficiently small, one has H p r\ > Co > where \rj\ < 2Sq, 
oj 1 / 2 < 2(e i5o) 1 ^ 2 - So H p 4> > on supp a, provided that \H p to\ < there. 

But being a sum of squares of functions with non-zero differentials, H p cu vanishes 
at u> — and satisfies \H p uj\ < Co; 1 / 2 . Due to ( |6.5[ ), we deduce that \H p u>\ < 
2C{e5) 1 / 2 . So \H p uj\ < a f holds if c -fe > 2C*(ec5) 1 / 2 , i.e. if e > C'S for some 
constant C > independent of e, Note that this constraint on e, i.e. that it 
cannot be too small, gives very rough localization: the width of the parabola at 
7] = — S is roughly it is very wide, and in particular insufficient to 

prove the propagation of singularities along the bicharacteristics. The reason is 
simple: our localizing function lu has no relation to H p , so we cannot expect a more 
precise estimate. Note, however, that the estimate is still non-trivial! Indeed, it 
shows that singularities propagate in the sense that £ cannot be an isolated point 
of WFsc(w). (We required e € (0,(5] beforehand, but in fact we could have dealt 
with e < //J, even if /i > 1, if we localized slightly differently.) 

We need to adapt to to H p to get a better estimate. If we linearize H p as above, 
and take uj = q 2 + . . . + ^In-ii then H p cu = and any e > works. Thus, in this 
case, we can prove propagation of singularities much like by the previous, simpler, 
construction. 

However, we do not need such a strong relationship to H p . Suppose instead that 
we merely get cu 'right' at £, in the sense that u> = J2 a j an( i Hp a j(£.) = 0- Then 
\H p <Tj\ < CoC^ 1 / 2 + |?7|), so |iJ p ^| < Cuj 1 ' 2 ^ 1 ' 2 + \t 1 \). Using , we deduce that 
\H p uj\ < fe provided that fe > C'ieSf^S, i.e. that e > C'8 3 for some constant 
C independent of e, S. Now the size of the parabola at r\ — —S is roughly uj 1 / 2 ~ S 2 , 
i.e. we have localized along a single direction, namely the direction of H p at £. By 
a relatively simple argument, one can piece together such estimates (i.e. where the 
direction is correct 'to first order') and deduce the propagation of singularities. We 
emphasize that the lower bound for e is natural. Indeed, with qj as above, we may 
take Uj e.g. to be <x, = qj + q\, j > 2. The bicharacteristics are qj — constant, 
but we are localizing near Oj — constant, and at rj = —5 these differ by S 2 . So 
any localization better than uj 1 / 2 ~ S 2 would in fact contradict the propagation of 
singularities! 

The microlocal positive commutator estimates in many-body scattering arise by 
this method. In particular, one can take V — which is the radial component 
of the momentum. The function to needs to be 7r-invariant, so if £ £ sc T*X a , 
it involves functions on sc T*X a as well as Z a = fjr and 77. The only additional 
argument needed is to show that the commutator is indeed positive, which has to 
be understood in an operator sense. Thus, the key point is that the commutator of 
H — A with a quantization B of rj is positive, modulo compact terms, when localized 
at A in the spectrum of H and microlocalized away from the radial set 7?.(A) . Note 
that, as usual, there is nothing to prove at 7?.(A), since each point in it is the image 
of a maximally extended generalized broken bicharacteristic. 

This positivity can be proved directly by showing that the indicial operators 
of the commutator are positive away from 1Z(\), which follows from an iterative 
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Figure 9. Bicharacteristics and suppa. The labels from Figure [s] 
have been removed to make the picture less cluttered. The straight 
horizontal lines are the a = constant lines, while the nearby 
parabolae are the bicharacteristics. 



argument. However, it also reduces to the Mourre estimate, involving the generator 
of dilations A — \(-D z + D z ■ z), which has principal symbol at £ • z. The Mourre 
estimate states the following. Suppose that A ^ A and e > 0. Then there is a 
8 > such that if <fi G (R) is supported in (A — 8, A + 8) then there exists 
K' G *~ 0O,1 (X,C) such that 

(6.6) <t>{H)i[A, H]<j>(H)) > 2(d(X) ~ e)<P{H) 2 + K>, 
where 

d(A) = inf{A - A' : A' < A, A' e A} > 
is the distance of A to the next threshold below it if A > inf A, and d(X) an arbitrary 
positive number if A < d(X). Since d(X) > if A ^ A, (6J3) is indeed a positive 
commutator estimate, which does not even have a 'negative' term E, unlike (6.1). 
The Mourre estimate, originating in ]47| , has been well understood since the work 
of Perry, Sigal and Simon [^3| and Froese and Herbst p6| . Here I just outline 
the argument in the simplest case, namely if no proper subsystem has any L 2 - 
eigenvalues. 

In this simplest case, the argument of |l6| (see also S and |37j) proceeds as fol- 
lows. In order to prove ( |6.6[ ), one only needs to show that for all b, the corresponding 
indicial operators satisfy the corresponding inequality, i.e. that 

(6.7) 4>(H b )i\AJI} b (t>(H b ) > 2(d(A) - e)dp{H b ) 2 . 

(This means that the operators on the two sides, which are families of operators on 
X b , depending on {y bl Q b ) G sc T^ b X, satisfy the inequality for all (yb, (b) € sc T^ b X.) 
It is convenient to assume that <f> is identically 1 near A; if ( |6.7D holds for such <fi, 
it holds for any 0o with slightly smaller support, as follows by multiplication by 
4>vi{Hb) from the left and right. 

Note that for b = the estimate certainly holds: it comes from the Poisson 
bracket formula in the scattering calculus, or from a direct computation yielding 

i[A, H] Q — 2A go . Hence, if the the localizing factor 0(i?o) = 0(IC| 2 ) ^ s supported 
in (A — 8, A + 8) and A > 0, then (6.7) holds even with d(X) — e replaced with A — 8. 
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Note that A > d(X), i f A > 0, since is a threshold of H. On the other hand, if 
A < 0, both sides of ( |6.7[ ) vanish for <f> supported near A, so the inequality holds 
trivially. 

In general, we may assume inductively that at all clusters c with C c C C b , i.e. 



X b C X c , (6/7) has been proved with <fr replaced by a cutoff cf> and e replaced by 
e', i.e. we may assume that for all e' > there exists 5' > such that for all c with 
C c C C b , and for all <j> € Cf (M; [0, 1]) supported in (A - 5', A + 5'), 

(6.8) 0{H c )i\A\H}J{H c ) > 2(d(A) - e')4>{H c ) 2 . 

But these are exactly the indicial operators of 0(if;,)i[^4, H] b (f>(H b ), so, as discussed 



in |6q, Proposition 8.2], (6/7) implies that 

(6.9) 0(^)ipC^6^(^6) > 2(d(A) - e')4>(H b ) 2 + K b , K b E % c °°' 1 (X b ,C b ). 

This implication relies on a square root construction in the many-body calculus, 
which is particularly simple in this case. 



Now, we first multiply (6.9) through by <f>{H) from both the left and the right. 
Recall that we use coordinates (z b ,z b ) on X b © X b and (Cb,C b ) are the dual co- 
ordinates. We remark that H b — \(^ b \ 2 + H b , so if A — \C, b \ 2 is not an eigenvalue 
of H b , then as supp0 — > {A}, <f)(H b + \( b \ 2 ) — > strongly, so as K b is compact, 
(f)(H b + \( b \ 2 )K b — > in norm; in particular it can be made to have norm smaller 
than e' — e > 0. After mu ltiplication from both sides by (f>i(H b ), with <f>\ having 



even smaller support, (6/7) follows (with <f>\ in place of 4>), with the size of supp^i 
a priori depending on However, i<fii(H b )[A, H] h 4>i{H b ) is continuous in ( b with 



values in bounded operators on L 2 (X b ), so if (6.7) holds at one value of ( b , then it 
holds nearby. Moreover, for large \( b \ both sides vanish as H b = H b + \( b \ 2 , with 
H b bounded below, so the estimate is in fact uniform if we slightly increase e > 0. 

In general, the proof requires to treat the range of E, the spectral projection of 
H b to {A}, separately. Roughly, the positivity estimate on the range of E comes 
from the virial theorem, iE[z b D z b, H b ]E = 0, which is formally clear, and is easy 
to prove. Thus, 

iE[A, H b ]E = iE[z b D z b , H b ]E + iE[z b D Zb , A Xb ]E - i[z b D Zb , A Xb ]E, 

and the commutator i(f>(H b )[z b D Zb , Ax b ]4>(H b ) is easily computed to be positive. 
Of course, there are also cross-terms that need to be considered, but they can be 



estimated by Cauchy-Schwartz estimates, see [|16| or |67 . 

I refer to [ |66| and |]65| for the detailed arguments proving propagation of singu- 
larities in the many-body setting, and to |38], Appendix] for weaker estimates with 
simplified proofs. 

7. Asymptotic completeness 

Asymptotic completeness (AC) is an L 2 -based statement describing the long- 
term behavior of solutions of the Schrodinger equation. In the short-range setting 
it says that for any <f> € L 2 (Xq) in the range of Id— E pp , E pp being the projection 
onto the bound states of H (i.e. onto its L 2 -eigenfunctions), there exist cf> a € L 2 (X a ) 
such that 

||e" iH V -^2e-* H '\4> a ® ip a )\\ -> as t -> +oo. 
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In the long-range setting, e~ tHat must be somewhat modified. After the ground- 
breaking work of Enss JlJ, [To), AC was first proved by Sigal and Soffer |55| in 
the short-range setting (see Graf's paper pof f or a different proof), and later by 
Derezihski and also by Sigal and Soffer |56|, [57| , in the long-range setting. In 
the short range case the main ingredient is equivalent to certain estimates of the 
resolvent at the radial sets in a sense that I now describe. In the long-range setting, 
as already in two-body scattering, additional constructive steps are needed. The 
estimates, in a different language, appeared first in the work of Sigal and Soffer 
(55). I hope that the following discussion makes it clearer how they relate to the 
propagation of singularities. 

While asymptotic completeness gives a complete long-term i 2 -description of 
solutions of the Schrddinger equation, the question remains whether an analogous 
description exists on other spaces, such as weighted L 2 -spaces. For example, if </> is 
Schwartz, are the <f> a Schwartz? Or dually, starting with a tempered distribution 
</>, are there tempered distributions a such that the convergence holds, as t — > +00, 
in a suitable sense? A different point of view is the parameterization of generalized 
eigenf unctions of H using the Poisson operators P Qj +(A), and the analogues of these 
questions can be asked there as well. The answer is affirmative in the two-body 
setting (even in the geometric setting, see jy], However, as indicated by the 

related issue of the mapping properties of the scattering matrices, discussed at the 
end of this section, it is unlikely that the same holds in the many-body setting. One 
can then ask weaker question, e.g. whether it holds in weighted spaces L 2 ., s near 
0. Or, one may ask whether one can give a precise description of the map tf> 1— ► <f> a 
e.g. as some sort of Fourier integral operator. 

As a starting point of relating the propagation of singularities to AC, we note 
that the propagation of singularities is proved by showing its 'relative' versions, i.e. 
that for any I, WF^ (u) is also a union of maximally extended generalized broken 
bicharacteristics. When considering the resolvent, first recall that for / G C°°(X), 
R(X + iO)f G H°°' 1 for all I < -1/2, so we only need to find WF*J(R(\ + i0)f) for 



I > -1/2. Theorem |3j is also valid for WF*J, i.e. the following holds 



Theorem 7.1. If X <£ A, then for f G 5(R"), / > -1/2, WF sc (i?(A + i0)f) is a 
subset of the image of 1Z+ (A) under the forward generalized broken bicharacteristic 
relation. 

This result allows u — R(X + i0)f not to lie in H*-^ 1 ' 2 on the image of TZ+(X) 
under the forward generalized broken bicharacteristic relation. This is a small set, 
but it is important to know whether WF*^ (u) may indeed intersect the forward 
image of 1Z+ (A). Of course, we cannot expect an improvement at 1Z+ (A), as shown 
already by the example of the free Euclidean Laplacian. The crucial improvement 
is the following estimate, due to Sigal and Soffer j55j. 

Theorem 7.2. If X ^ A, then for f G S(R n ), WF*-" 1/2 (i?(A + i0)f) C K+(X). 

Remark. This theorem also has a time-dependent analogue. If it is a solution of 
the Schrodinger equation (D t + H)u = with u\ t= o G S(Xo) then on the one hand 
u G H°° l (X x K) for I < -1/2, on the other hand WF*-~ 1/2 (m) C K. 

In fact, this theorem can be improved along the lines of the distributional state- 



ment in Theorem B.E 
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Corollary 7.3. Suppose that A $ A, f e H*' 1/2 and WF*;V 2+£ (/) n Tl-(X) = 
for some e > 0. Then R(X + i0)/ = lim t -o #(A + it)/ eiiste in ff* -i/2-e' (x), 
e' > 0, and WF*^ 1/2 (i?(A + i0)/) C U + {X). 



Theorem 7.2 can be proved rather simply. The main issue is how to obtain a pos- 
itive commutator at the radial point. Away from the radial sets arbitrary weights 
can be accommodated by suitable construction, as pointed out in the previous sec- 
tion. At radial points only the weights can give positive commutators. Now, one 
has to use weights x^ 21 ^ 1 to obtain estimates for WF*^, and these weights will give 
a commutator whose sign depends on that of —2/ — 1, hence on whether / > —1/2, 



I < —1/2 or I — —1/2. It turns out that the sign is correct for (6.3) to be of 
use if I < —1/2; this yields the limiting absorption principle. The sign is wrong if 
I > —1/2, so no results can be expected then. In the borderline case I = —1/2, 
the weight vanishes. The way to obtain a positive commutator is thus to consider 
operators A which are microlocally (a multiple of) the identity near 1Z + (X). The 
commutator then vanishes microlocally near 1Z + (A), which is reasonable since no 
estimate on WF*; _1 ' 2 (ti) can be expected there. 



It is then straightfoward to construct A so that (|6.3| ) can be used to prove 
Theorem 7.2 . Indeed, it suffices to show that on WF* ( ;~ 1 ^ 2 (u), V — must satisfy 



A — rj 2 e A, for then the full statement of the theorem follows by the propagation of 
singularities for WF* c _1 ' 2 (m). So we proceed to prove this simpler result, namely 
that if A - f\ 2 £ A then for any point £, 77(f) = 77 implies f </ WF*,: -1 / 2 ^). 

To do so, we let a — x(v) where \ € C£°(R), X > 0> is chosen so that x = 1 
on [0,77 — S] for some S > 0, x' < on (0,oo), x'(v) < 0, and t 6 suppx' implies 
that A - t 2 A. This can be arranged as A is closed. We can further make sure 
that \J — x' is C°° on (0, 00). Then the positive commutator methods outlined show 
the commutator of A, a quantization of a, with H — A is positive, in the region 
77 > 0, yielding the estimate that proves the theorem. We remark that partial 
microlocalization, using functions of 77, hes been used extensively in many-body 
scattering, especially by Gerard, Isozaki and Skibsted [|l8| and Wang to 
obtain partially microlocal statements such as radiation conditions and uniqueness 
statements |3^, |3l]], and indeed to prove the smoothness of 2-cluster to 2-cluster 
scattering matrices Q . 



It turns out that there is an even simpler way of proving Theorem 7.2, or indeed 
a stronger statement, which is due to Yafaev |Q. His estimate states that in a 
neighborhood of C atteg , where we write y a for the coordinates z a /\z a \ along C a . lcg , 
xD Va R(\ + i0)f is in H*'^ X / 2 {X). Since the principal symbol of xD Va is invertible 



on (T*X a) re g H Char(A)) \ 7£(A), this result implies Theorem 7.2. Yafaev's proof 
relies on a simple and explicit commutator calculation, which allows one to deal 
with various error terms that one may, a priori, expect. However, exactly because 
of its explicit nature, it is presumably hard to generalize to more geometric settings, 
while the argument we sketched does not face this difficulty. 



As discussed by Yafaev [73 in the usual time-dependent version, short-range 



asymptotic clustering, hence asymptotic completeness, are relatively easy conse- 



quences of Corollary 7.3, and we refer to |73J and |S| for more details. However, it 
is worth pointing out that the reason why Coulomb-type potentials (i.e. those in 
S^ 1 ) are not 'short-range' is that the Hamilton vector field in some subsystem van- 
ishes at radial points. This degeneracy makes even the subprincipal term important 
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in describing the precise behavior of generalized eigenfunctions microlocally near 
this point. 

Before turning to scattering theory on symmetric spaces, we note the implications 



of Theorem 7.2 for the scattering matrices. Previously, S a p(X) was only defined as a 
map S a p(X) ■ C^°(S a . Icg ) — > C~°°{Sb,xeg)- Indeed, part of the broken bicharacteristic 
relation connects 1Z+(X) with its image, and this can a priori give a singularity in 
the kernel of S a p(X) of the kind that does not even allow one to conclude that 
S a p(X) : C^°(Sa :IC g) — > C°°(Sb,r C g)- The pairing formula, (4.2), combined with 



Theorem 7.2, show that in fact 

(7.1) S a p(X) : L 2 (S a ) L 2 (S b ). 

It is an interesting question whether this can be improved if we restrict S a p(X) to 
(S a , rcg)- Namely, except in special cases such as TV-clusters and two-clusters, 
the best known result is the trivial consequence of ( [T.l]) : 

S a p(X) : C^°(S a ^ cg ) — > L 2 (St)- 

(In the case of iV-clusters and 2-clusters, the geometry of generalized broken bichar- 
acteristics gives S a p(X) : C^°(S a ^ cg ) — > C°° (Sb, le g) ■) The putative improvement 



would have to be connected to an improvement of Theorem 7.2, namely to the 
existence of some I > -1/2 such that WF*^(R(X + iO)f) C 1Z+ (A). It would also 
be connected to a better understanding of R(X ± iO) at the thresholds, in which 
direction Wang's paper Q is the only one I am aware of. 



8. Scattering on higher rank symmetric spaces 

In this section I discuss SL(N, R)/ SO(N, R), indeed, mostly SL(3, R)/SO(3, R). 
The books p5| , p4[ and are good general references. N — 2 yields the hyperbolic 
plane H 2 , which is a rank one symmetric space on which many aspects of analysis, 
such as the asymptotic behavior of the resolvent kernel and the analytic continu- 
ation of the resolvent are well understood. Indeed, these have been described on 
asymptotically hyperbolic spaces by Mazzeo and Melrose [|7| and Perry H|| . 

Higher rank symmetric spaces, such as SL(N)/ SO(iV), N > 3, are much less un- 
derstood. For example, using results of Harish-Chandra, and Trombi and Varadara- 
jan (see pf|), Anker and Ji only recently obtained the leading order behavior of the 
Green's function ^, H, 0). Also, while spherical functions, which are most analogous 
to partial plane-partial spherical waves in the Euclidean setting, have been analyzed 
by Harish-Chandra, Trombi and Varadarajan, and in particular their analytic con- 
tinuation is understood, the same cannot be said about the Green's function. The 
analysis of spherical functions relies on perturbation series expansions, much like 
in the proof of the Cauchy-Kovalevskaya theorem, and it does not work well at the 
walls of the Weyl chambers. Here I only illustrate some recent joint results with 
Rafe Mazzeo [^0^ that illuminate the connections with many-body scattering, 
and in particular give rather direct results for the resolvent. 

First I describe the space SL(3) / SO(3). The polar decomposition states that any 
C G SL(3) can be written uniquely as C = VR, V = (CC*) 1 / 2 is positive definite 
and has determinant 1, R G SO(3). Thus, SL(3)/SO(3) can be identified with 
the set M of positive definite matrices of determinant 1; this is a five-dimensional 
real analytic manifold. The Killing form provides a Riemannian metric g. The 
associated Laplacian A = A g gives a self-adjoint unbounded operator on L 2 (M, dg), 
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with spectrum [Aq,+oo), Ao = |. Let R(X) = (A — A) 1 be the resolvent of A g , 
A g [A ,+oo). 

Fix a point o £ M, which we may as well assume is the image of the identity 
matrix / in the identification above. The stabilizer subgroup K Q (in the natural 
SL(3) action on M) is isomorphic to SO(3). The Green function G Q (A) with pole 
at o and at eigenvalue A is, by definition R(X)5 . It is standard that G a lies in 
the space of i^ -invariant distributions on M. It is thus natural to study A on 
-Ko-invariant functions. 

Perhaps the most interesting property is the analytic continuation of the resol- 
vent, which I state before indicating how it, and other results, relate to many-body 
scattering. 

Fix the branch of the square root function J on C \ [0, +oo) which has negative 
imaginary part when w £ C \ [0, +oo). Let S denote that part of the Riemann 
surface for A i— ► yA — Ao where we continue from A — Ao ^ [0, +oo) and allow 
arg(A — Ao) to change by any amount less than tt. In other words, starting in the 
region Im y/X — Ao < 0, we continue across either of the rays where Im y/X — Ao = 
and Re V A — Ao > 0, respectively < 0, allowing the argument of VA — Ao to 
change by any amount less than tt/2 (so that only the positive imaginary axis is 
not reached). 

Theorem 8.1. With all notation as above, the Green function G (X) continues 
meromorphically to S as a distribution. Similarly, as an operator between appro- 
priate spaces of K -invariant functions, the resolvent R(X) itself has a meromorphic 
continuation in this region, with all poles of finite rank. 

Having stated the theorem, I indicate how it relates to many-body scattering. To 
do so, fix the point o - we may as well take it to be the identity matrix /. Now, M is 
a perfectly nice real analytic manifold and A is an elliptic operator on it in the usual 
sense, so the only question is its behavior at infinity. In order to describe this, we 
remark that any matrix A £ M can be diagonalized, i.e. written as A — OAO', with 
O € SO (3) and A diagonal, with entries given by the eigenvalues of A. If a is the set 
of diagonal matrices of trace 0, then A £ exp(a). If all eigenvalues of A are distinct, 
then A is determined except for the ordering of the eigenvalues, and there are only 
finitely many possibilities for O as well. However, at the walls, which are defined 
to be the places where two eigenvalues coincide, there is much more indeterminacy. 
For example, if two eigenvalues coincide, only their joint eigenspace is well-defined. 
Correspondingly, we may replace O by O'O for any O' £ SO (3) preserving the 
eigenspace decomposition and still obtain the desired diagonalization. 

This is closely reflected in the structure of the Laplacian at infinity. In fact, it 
turns out that on SO(3)-invariant functions, A is essentially a three-body Hamil- 
tonian on a with first order interactions and with collision 'planes' given by the 
walls (they are lines), see e.g. |25|, Chapter II, Proposition 3.9]. So rather than 
particles, eigenvalues scatter in this case! Consequently, many-body results can be 
adapted to this setting. 

We indicate how this is done. The most succint way of describing the geometry of 
M at infinity is to compactify it to a manifold M with codimension two corners. It 
has two boundary hypersurfaces, ffj and H\ which are perhaps easiest to describe 
in terms of a natural system of local coordinates derived from the matrix represen- 
tation of elements in M. As above, we write A £ M as A — OAO*, with O £ SO(3) 
and A diagonal. The ordering of the diagonal entries of A is undetermined, but in 
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the region where no two of them arc equal, we denote them as < Ai < A2 < A3 
(but recall also that A1A2A3 = 1). In this region the ratios 

Ai A2 

are independent functions, and near the submanifold cxp(a) in M we can complete 
them to a full coordinate system by adding the above-diagonal entries C12, C13, C23 
in the skew-symmetric matrix T — logO. On M we have /i, v > 0, and locally the 
compactification consists of replacing (fi, v) S (0, 1) x (0, 1) by (fi, v) S [0, 1) x [0, 1). 
Then H} = {/i = 0} and H$ — {v — 0}, and this coordinate system gives the C°° 
structure near the corner ify PI if 




H t na 



Figure 10. The closure of a, or rather exp(a), in the compacti- 
fication M of M. The lines in the interior are the Weyl chamber 
walls, playing the role of collision planes in many-body scattering. 
The side faces H* n a and H$ fl a correspond to the front faces on 
Figure |. The main face on Figure |5| would only show up if we 
did a logarithmic blow-up of all boundary hypersurfaces of M and 
then blew up the corner. 

On the other hand, in a neighborhood of the interior of H*, for example, we 
obtain the compactification and its C°° structure as follows. Write the eigenvalues 
of A <E M, i.e. the diagonal entries of A in the decomposition for A above, as Ai, 
A2 and A3 . Suppose that A lies in a small neighbourhood U where 

c<Ai/A 2 <c _1 , A 3 > 1/c, 

for some fixed c € (0, 1). Recall also that A3 = I/A1A2. These inequalities imply 
that Ai = (Ai/A 2 ) 1/2 A3~ 1/2 < 1 and A 2 = (A 2 /Ai) 1 / 2 A 3 _1/2 < 1, and A 3 > 1 in U. 
Hence there is a well-defined decomposition R 3 = En © E3 for any A 6 U, where 
E12 is the sum of the first two eigenspaces and E3 is the eigenspace corresponding 
to A3, regardless of whether or not Ai and A 2 coincide. We could write equivalently 
A = OCO 1 , where C is block-diagonal, preserving the splitting M 2 © R of M. 3 . The 
ambiguity in this factorization is that C can be conjugated by an element of 0(2) 
(acting in the upper left block), and 0(2) can be included in the top left corner 
of SO(3) (the bottom right entry being set equal to ±1 appropriately). Let C 
denote the upper- left block of C; the bottom right entry of C is just A3, and so 
AsdetC" = 1. In other words, C — A 3 C" where C is positive definite and 
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symmetric with determinant 1, hence represents an element of SL(2)/ S0(2) = H 2 . 
Hence for an appropriate neighbourhood V of / in SL(2)/ S0(2), the neighbourhood 
U is identified with (V x SO(3))/0(2) x (0,c 3/2 ), where the variable on the last 
factor is s — A 3 3 ^ 2 . The compactification then simply replaces (0, c 3 / 2 ) s by [0, c 3 / 2 ). 
Note that although the action of 0(2) on V has a fixed point (namely I), its action 
on SO(3), and hence on the product, is free. The neighbourhood V can be chosen 
larger when A3 is larger, and the limiting 'cross-section' A3 = const has the form 
(H 2 x SO(3))/0(2). This space is a fibre bundle over SO(3)/0(2) (= RP 2 ) with 
fibre H 2 . Notice that the Weyl chamber wall corresponds to the origin (i.e. the point 
fixed by the SO(2) action) in HI 2 . I refer to |I(J f° r a more thorough description of 
M. 

On each boundary hypersurface of M, it is now easy to describe model operators 
for A acting on SO(3)-invariant functions. For instance, at H$ this model can be 
considered as an operator Lj on R s x H 2 , acting on SO(2)-invariant functions. 
Explicitly, 

h = \(sD s ) 2 + i\{sD s ) + ^A m . 

This is tensor product type, so its resolvent can be obtained from an integral of the 
resolvents of j(sD s ) 2 + i^(sD s ) and gAjp. (Note that I am ignoring the weights 
of the L 2 spaces on which we are working, hence the appearance of the perhaps 
strange first order terms.) 

This framework allows one to develop the elliptic theory, for example to analyze 
(A — A) -1 for A € C \ [Ao,+oo). In particular, one can construct a parametrix 
for A on M that has a smoothing error. Since this error has no decay at infinity, 
it is not compact. However, the error can be improved by pasting together the 
resolvents of Lj and L", and applying the result to the error term to remove it 
modulo a decaying, hence compact, new error term. One of the consequences is 
then the description of the asymptotic behavior of the Green's function, see Eot]. 



The proof of Theorem 8.1 relies on the method of complex scaling. This is an 
extension of dilation analyticity, and was originally developed in the setting of 2- 
body scattering by Aguilar-Combes and generalized to the m any -body setting 
by Balslev-Combes || and later further generalized by Gerard fl7|] . We refer to 
(26) and |34[ Volume 4] for an exposition, and to the paper |38) of Sjostrand and 
Zworski for a slightly different point of view. 

The point of complex scaling is to rotate the essential spectrum of the operator 
being studied, in this case the Laplacian. To give the reader a rough idea how 
this works, consider the hyperbolic space H 2 = SL(2, R)/ SO(2, R), which may be 
identified with the set of two- by- two positive definite matrices A of determinant 1. 
In terms of geodesic normal coordinates (r, to) about = I, the Laplacian is given 

by 

A H 2 = D 2 r . - i(cothr)D r + (sinhr) -2 !) 2 ,. 

Now consider the diffeomorphism <& e : A i-» A w , w = e 6 , on H 2 , 6 e R. This 
corresponds to dilation along the geodesies through o, since these have the form 
■ s 1— > A cs , c > 0. Thus, in geodesic normal coordinates, <&e : (r,u>) 1— > (e e r, uj). 
$e defines a group of unitary operators on L 2 (H 2 ) via 

(Ugf)(A) = (det£>A*e) 1/a («5/)(^), J = detD A ^ 9 = w^^, w - e e . 

smnr 
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Now, for 9 real, consider the scaled Laplacian 

(A H2 ) e = U e A m ,U^ = J^&gAw&gJ' 1 ' 2 

= J 1/2 {w- 2 D 2 r - iw- 1 coth(wr)D r + (sinh^r))- 2 ^ ) j-1/2 

This is an operator on H 2 , with coefficients which extend analytically in the strip 
|Im0| < If. The square root is continued from the standard branch near w > 
0. (The singularity of the coefficients at r = is only an artifact of the polar 
coordinate representation.) Note that (A H 2) e and (Ajjp)^ are unitary equivalent if 
lm6 = lm9' because of the group properties of Ug. The scaled operator, (A$p)g, 
is not elliptic on all of H 2 when < |Im#| < § because for r large enough, 
w 2 sinh(wr)~ 2 can lie in R~. However, it is elliptic in some uniform neighbourhood 
of o in H 2 , and its radial part 

(A H 2) 0irad = J x l 2 {w- 2 D 2 r - iw- 1 coth(u,r)A-y~ 1/2 , 

which corresponds to its action on SO(2)-invariant functions, is elliptic on the entire 
half-line r > 0. The model operator for (Ajj2)g ra( j — A at infinity, 

e (t»-l)r/2 (u; -2 £ ,2 _ lw -l Dr _ x)e -{ W -l)r/2 

= e ( - w - 1)r/2 ((w~ 1 D r - l -) 2 - (A - I)) e -(— !W 2 , 

is also invertible on the model space at infinity, L 2 (R; e r dr), since this is equivalent 
to the invertibility of 

e-'\{ W -'D r \) 2 - (A - i))e— /» = W ~ 2 D 2 -(X-\) 

on L 2 (R;dr). Thus, a parametrix with compact remainder can be constructed for 
(A H 2) e rad , and this show that its essential spectrum lies in j + e~ 2lIm0 [O, +00). 
Hence ((A H 2)g jra d — A) -1 is meromorphic outside this set. In fact, it is well known 
that there are no poles in this entire strip (although there are an infinite number 
on arg V A — Ao = 7r/2). 

Combining this with some more standard technical facts, we are in a position to 
apply the theory of Aguilar-Balslev-Combes to prove that ((An2) ra d — an d 
hence (A H 2 — A) -1 , has an analytic continuation in A across ( \, +00). This is done 
by noting that for SO(2)-invariant functions /, g £ L 2 {M 2 ) and 9 eR, 

(/, ((A H 2) rad - X)~ 1 g) = (Ugf, ((A H 2) , rad - A)" 1 ^) 

by the unitarity of Ug. Now if /, g lie in a smaller (dense) class of functions such that 
Ugf and Ugg continue analytically from figR, then the meromorphic continuation 
in A of the right hand side is obtained by first making 9 complex with imaginary 
part of the appropriate sign, and then allowing A to cross the continuous spectrum 
of A H 2 without encountering the essential spectrum of (A H 2)e rad . Hence the left 
hand side continues meromorphically also. With some additional care, one can even 
allow g to be the delta distribution at o, yielding the meromorphic continuation of 
the Green's function. 

The argument on the higher rank symmetric space M = SL(3)/ SO(3) is similar. 
We still use the same scaling $g : A A w , w = e e , on M, S £ M. Again, the 
first concern is that, allowing 9 to become complex, the scaled operator Ag is not 
elliptic. However, it is elliptic near o — Id, and the scaled models for it near the 
walls, such as {L$)g, remain elliptic at the walls. After all, for the latter, this is 
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just the ellipticity of (Anp)g near the origin, which we have already observed. This 
again allows the elliptic parametrix construction to proceed, supplying the results 
we needed in order to reach the framework of complex scaling. This in turn finishes 
the proof of Theorem 
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